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1 - INTRODUCTION 


It is widely understood that the use of mathematical modeling, analysis, and 
simulation (MMA&S), when used wisely, is of great value as a means to gain early 
insights into the nature of a proposed engineering system and to enhance the 
probability of anticipating problem areas early in the program. It is also now 
widely understood that MMA&S, when properly integrated into a system 
development program, can result in a system development cost which is orders of 
magnitude less expensive than the cost of development via the once traditional 
strategy of rushing through the preliminary design in an ad lib fashion to the early 
fabrication of the physical system, testing, changing design, rebuilding, retesting, 
etc. A combination of MMA&S, design, fabrication, and testing will minimize total 
system development cost. Systems control is an area of one of the earliest 
applications of MMA&S, possibly because of the subtle and often counter-intuitive 
nature of the dynamic stability problem. Results have been extraordinarily 
successful in many cases. 

The Controlled Ecological Life Support System (CELSS) concept presents a 
prospective application of this system development strategy. It is interdisciplinary 
in its scientific and technological span and, of necessity, contains dynamic processes 
of wide diversity of physical origin. Such a system would contain at least a crop 
growth subsystem, a personnel subsystem, a waste treatment subsystem (Figure 1), 
and a management, automation, and control subsystem. Systems characterized by 
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this level of complexity and diversity almost always have the potential for counter- 
intuitive behavior which, at best, results in unacceptable performance of the 
controlled system, and, at worst, results in self-destructive behavior when 
im properly chosen feedback control is imposed on the system. 



Figure 1. Schematic showing the primary functional subsystems of a CELSS. Data 
acquisition, data processing, and control subsystems are not shown. 

The Crop Growth Research Chamber (CGRC) project, an element of the 
CELSS program, can itself be profitably subjected to the MMA&S treatment. 
(As pointed out earlier, CELSS must contain a crop growth subsystem, but its 
required envelope of operating conditions is much smaller than that of a 
CGRC.) The CGRC is intended to be a precision research apparatus to support 
research in crop growth science and technology and must provide a precisely 
controlled environment for this purpose. The performance requirements include 
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providing satisfactory (prescribed) time histories of fluid flows, concentrations 
of gases and liquids, temperatures, and pressures at designated points within the 
CGRC system, a scenario which is referred to as the multi-input, multi-output 
tracking problem in the field of systems control. Even a subsystem of a CGRC, 
such as the shoot side environmental control system, is sufficiently complex to 
present a challenging but clearly tractable and almost ideal subject for MMA&S 
attention. 

The development of a system to control the environmental conditions 
within the shoot side chamber of the CGRC involves (1) the selection of 
component devices (i.e. actuators) which will affect the environmental variables 
that are to be controlled, (2) the selection of component devices which will 
provide suitable measurements of the environmental variables (i.e. 
measurement instruments), and (3) the development of control algorithms 
which utilize the outputs of the measurement instruments to provide the 
necessary combinations and time sequencing of inputs to the actuators so that 
the required performance characteristics of the environmental variables within 
the CGRC are met. 

The process of control system synthesis almost always involves a 
familiarization process, and the familiarization process almost always involves a 
process of starting with great but meaningful simplifications and moving in 
some rational sequence of stages to a stage in which sufficient complexity and 
detail is addressed. The conscious pursuit of this strategy and what it implies is 
apparently not generally intuitive. 

The purpose of this report is to present the results of the authors’ initial 
examination of the control problem related to the provision of the required 
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environment of the plant shoots of a hypothetical CGRC configuration. This 
configuration was thought to be consistent with the requirements of a CGRC 
but not expected to be the final configuration. Figure 2 is a minimally 
annotated schematic representation of the object of study, and is provided to 
help the reader gain insights into the relative positions of the major components 
and their sequence in the air flowstream. 


1.1 ORGANIZATION OF THE REPORT 

In the sequel, we describe the crop shoot side environmental delivery 
system (CSSEDS), define the specific purpose of the model, characterize the 
model in terms of the processes of which it is composed and how they interact, 
present the natural equations in both linear and nonlinear forms, extract the 
state variable representation of the model, present a canonical form of the state 
variable form, and present the results of a cursory study of the problem of 
choosing feedback control gains so that one of the unstable modes of the system 
is replaced by a stable mode. 

2 - DESCRIPTION OF THE SHOOT ENVIRONMENT DELIVERY 

SYSTEM 

Figure 3 is also a representation of the CSSEDS but is configured to be 
very useful to the modeling process and displays much more detail than does 
Figure 2. In Figure 3, large circled numbers indicate points (nodes) at which it 
is anticipated that the potential variables are important to the control and 
performance problem, and small numbers beside arrows with 
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Figure 2. Diagram of the crop shoot side environment delivery system. Relative position 
of primary functional subsystems and fluid flow paths are empha sized. 

one-sided heads indicate flow paths for which it is anticipated that the flow 
variables will have similar importance. The function of the components of the 
CSSEDS will be explained by tracing the airflow clockwise around the circuit 
beginning with its entrance into the upper chamber. The plants receive radiant 
energy from a light source above (and external to) the chamber after it emerges 
from an electromagnetic radiation biter which removes most of the long-wave (IR) 
radiation from the spectrum. Water and nutrients are supplied to the roots by 
means of a nutrient delivery system that is materially isolated from the shoot side 
environment. The flow path for the air entering the chamber is assumed to be 
designed to assure that uniform conditions exist in the shoot side c ham ber 

Air enters near the top of the upper chamber and is thoroughly mixed with 
the air in the upper portion of the chamber. As the air moves from the upper 
chamber inlet to its exit, gaseous exchange of carbon dioxide, water vapor and 
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oxygen occurs between the crop canopy and the atmosphere. The resulting mixture 
exits the upper chamber via a flow path provided between the chamber walls and 
the plant support surface into the lower chamber and then into the ductmg via an 
entrance located in the bottom. The flow is then directed through a pre-filter and a 
HEP A filter to remove particulates from the air as it leaves the chamber. Upon 
exiting the filter, the air sustains a pressure loss imposed by a valve. By means of 
a valve and fan, a portion of the air flow is diverted into a gas separator which 
removes excess oxygen and/or carbon dioxide as required. The required air 
movement around the flow path is accomplished by employing an air mover (e.g., a 
centrifugal pump or a blower), and at the discharge of the air mover, makeup gases 
are injected into the flow stream to maintain the required atmospheric 
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composition. The water which is transpired by the plants is removed from the 
system airflow stream in the form of condensate by diverting a portion of the flow 
through a condensing heat exchanger. Two valves, one in the flow path through the 
condensing heat exchanger and one in the condensing heat exchanger bypass 
regulate the mass flow ratio of the paths. The flow is recombined and enters the 
cooling/heating subsystem. A portion of the flow can be diverted through either a 
heater or a cooling heat exchanger. Valves control the flow distribution among the 
three possible flow paths. Again, the flows are combined and the mixture is ducted 
to the chamber inlet, completing the traverse of the circuit. 

3 - PERFORMANCE AND PHYSICAL CONFIGURATION REQUIREMENTS 

It was established by a science advisory group that a CGRC which could be 
useful for their research would of necessity have certain performance capabilities 
and would have certain physical characteristics. The group limited itself to 
describing these capabilities and characteristics and left the details of how they 
might be obtained to the designers. 

(1) . Crop growth platform area - 2.0 m 2 

(2) . Minimum upper crop shoot volume - 2.0 m 3 . 

(3) . With the lights on and no plants in the CGRC, the temperature rise of the air 

exiting the upper chamber must not be greater than 2 C. 

(4) . The CGRC must be able to provide a quantum light intensity of as much as 

2000 pmoles/(m 2 -sec). 

(5) . The airflow speed must not exceed 2 m/s in the shoot upper chamber 
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(7) . It is not permissible for the pressure within the shoot and root environments to 

be less than 0.5 mm Hg higher than the CGRC surroundings. 

(8) . The chamber walls must be kept at least 2 C higher than the air within the 

chamber. 

(9) . Water must not precipitate from the air anywhere other than in the water 

removal subsystem. 


4 - PURPOSE OF THE MODEL 

We consider the scenario in which the growth chamber contains a crop of 
plants forming a closed canopy and the system goal is to maintain temperature, 
relative humidity, pressure, carbon dioxide concentration, oxygen concentration, 
and m ean air velocity in the chamber within set tolerances about constant 
operating points. The formal purpose of the model is to support the exploration of 
the generic stability properties of the system in the neighborhood of an operating 
condition of specific interest to members of a science advisory group, but it also 
played an important role in the process of familiarization. Of course, the scientists 
who expect to do research with this type of chamber are concerned about a great 
deal more than simple stability, but stability is a logical starting point for 
performance studies. The details of the operating condition for which the study was 
done are presented in Sections 3, 7.4, and 7.5. 
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5 - THE MATHEMATICAL MODEL OF THE CROP SHOOT ENVIRONMENT 

DELIVERY SYSTEM 

The philosophy and procedure with which the mathematical model of the 
CGRC system was developed is documented in [1] (see Appendix Al). Here, we 
present a summaxy of the procedure at the level of detail necessary for purposes of 
continuity of thought. 

The initial stage of the modeling process, the conceptualization of (1) which 
processes within the system are to be addressed, (2) how they interact with one 
another, (3) how the system interacts with the surroundings, and (4) how each 
process is to be characterized mathematically /numerically is a purely abstract 
process. To make each of these decisions tangibly available for the modeler’s 
utility, as well as the equally important responsibility of communicating the details 
of the model to others, preparing an adequately annotated list suffices effectively 
for the first purpose, constructing a “modeling system” diagram is very useful for 
the second and third purposes, and listing the “natural equations” of the system is 
very effective for the fourth. 

The behavior of any system is the resultant of the nature of the processes 
contained by the system, their interaction, the conditions at the time the system is 
set free, and the interaction of the the system with its surroundings. The use of 
mathematical modeling to predict the behavior of the system in the context of 
some question or questions requires (1) the recognition that these processes are in 
action, that they interact with one another and the details of those interactions, 
and that they interact with the surroundings and the details of those interactions, 
and (2) committing this recognition to a prescription of mathematical operations 
and parameter values. The intellectual process of accomplishing these actions is 
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called mathematical modeling. 


In the course of generating a system representation, one synthesizes an 
idealized conceptual system (the “modeling system”) of idealized processes 
(“modeling processes”) which mimic the actual system’s physical processes to an 
adequate degree. In this procedure, it is implicit that by contemplation of the 
actual system, one must correctly infer how the idealized processes must be 
arranged to interact so that the modeling system will emulate the actual system s 
behavior in the context of the problem being addressed to a degree of adequate 
accuracy. In Appendix A.l the reader can find a set of facts and thoughts about the 
CSSEDS which shape the contemplation. Logically, one first identifies the processes 
and then establishes how they interact. 

5.1- SUMMARY OF THE PHASES OF THE MODEL GENERATION 
PROCEDURE 

In the sequel, the description of the stages of the model generation 
procedure will be found; 

(1) List the physical processes which axe to be represented by model processes 

(2) Declare the interpretation of how the processes interact, sign conventions, 
point identities and path identities. This is done via process interaction 

diagrams Figures 4 and 5. 

(3) Write the primitive equations; the mathematical model of each process in 
terms of the point variables and path variables associated with that process 
in Figures 4 and 5, and the compatibility equations (those which express 
constraints among the point and path variables). In this phase, the linearized 
form of the equation is also developed. 


10 


(4) Solve the primitive equations for the nominal values of the variables at the 
experimental operating condition. 

(5) Write the linear primary equations. Take advantage of opportunities to 
eliminate variables ad lib using some of the equations. The object of this 
action is to reduce the number of equations which must be manipulated. 

(6) Substitute nominal values into linearized equations to obtain equation forms 
which have numerical coefficients. This form is characterized by the lineax 
equations: 


/4j77 = A 2 r) + A 3 <t + A 4 u + Acd 


(5.1.1) 


and 


Bj<t = B 2 rj + fi 3 u + B 4 d 


(5.1.2) 


The state variable format can then be determined, using, e.g., 'MATLAB 


x = Af \A 2 + A 3 B~ 1 B 2 )x + A{- \A 4 + A z B{ l B z )u + A~ ! (>1 5 + A 3 B~ l B 4 )d (5.1.3) 

and 

a = 1 (B 2 x+ B 3 u + B 4 d) (5.1.4) 


(generalized inverse is implied) which is of the simple form 


and 


x = J.x + Sa + C d 


(5.1.5) 


a — + 8u + ‘Jd 


(5.1.6) 


(7) Analyze the state variable representation for the purpose of familiarization 
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with its dynamical properties and its controllability properties. There is a 
great variety of choices of what might be done and it varies with the 
application and the investigator. In the sequel, we present the results of: 

1. the Jordan canonical transformation 

2. the Real Jordan transformation 

3. eigenvalue assessment 

4. controllability assessment, and 

5. an elementary feedback control study. 

5.2 - THE MODELING PROCESSES. 

It was anticipated that the following physical processes would have to be 
included in the modeling system in order for the representation to be useful: 

(1) . Fluid inertance in upper duct. 

(2) . Radiant energy flow from the light source into the upper shoot chamber 

(source of PAR). 

(3) . Molar storage of N 2 in upper chamber. 

(4) . Molar storage of 0 2 in the upper chamber. 

(5) . Molar storage of C0 2 in upper chamber. 

(6) . Molar storage of H 2 0 vapor in the upper chamber. 

(7) . 0 2 flux from crop into the upper chamber air. 

(7) . C0 2 flux from the upper chamber air into the crop. 

(8) . H 2 0 flux from the crop into the upper chamber air. 

(9) . Upper chamber gas state. 

(10) . Thermal energy of gas in upper chamber. 

(11) . Volume change of the upper chamber due to glove port activity. 
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(12) . Heat exchange between walls of upper chamber and upper chamber air. 

(13) . Heat exchange between crop canopy and upper chamber air. 

(14) . Flow energy loss; airflow past plant platform. 

(15) . Molar storage of N 2 in lower chamber. 



(16) . Molar storage of 0 2 in the lower chamber. 

(17) . Molar storage of C0 2 in lower chamber. 

(18) . Molar storage of H 2 0 in the lower chamber. 
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(20) . Flow energy loss; exit from lower chamber. 

(21) . Flow energy loss; filter and lower valve. 

(22) . Extraction of excess oxygen; lower duct section. 

(23) . Fluid inertance; lower duct section. 

(24) . Flow energy loss; lower duct section. 

(25) . Energy conversion by the pump. 

(26) . Motor and pump rotor rotational inertia. 

(27) . Control moment on the motor armature. 

(28) . Motor and blower friction. 
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(29) . Injection of C0 2 into the gas addition mixing section. 

(30) . Molax storage of N 2 in the gas addition mixing section. 

(31) . Molax storage of 0 2 in the gas addition mixing section. 

(32) . Molax storage of C0 2 in the gas addition mixing section. 

(33) . Molax storage of H 2 0 in the gas addition mixing section. 

(34) . Thermal energy storage in the gas addition mixing section. 

(35) . Flow energy loss; condensing heat exchanger bypass control valve. 

(36) . Flow energy loss; condensing heat exchanger air side flow control valve. 

(37) . Liquid water precipitation; condensing heat exchanger 

(38) . Thermal energy exchange; condensing heat exchanger. 

(39) . Flow energy loss; duct from condensing to heating/cooling system. 

(40) . Flow energy loss; air heater flow control valve. 

(41) . Flow energy loss; air cooler flow control valve. 

(42) . Flow energy loss; air heater/air cooler bypass flow control valve. 

(43) . Thermal energy exchange; air cooler heat exchanger. 

(44) . Flow energy loss; heater flow control valve. 

(45) . Volume change of the shoot upper chamber due to glove port activity. 

The entire CSSEDS is assumed to be adiabatic with respect to its 
surroundings. All ducts are assumed to be of the same cross-sectional area for flow. 

5.3 - THE MODELING SYSTEM 

In the coherent process of generating a system representation, one 
synthesizes an idealized conceptual system, the modeling system, of idealized 
processes. In the contemplation of the actual system, one must correctly infer 
which processes axe at work and how the idealized processes must be arranged to 
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interact so that the modeling system will emulate the behavior of the actual system 
in the context of the problem being addressed to a degree of adequate accuracy. 
Logically, one first perceives the processes that are in action and then how they 
interact. Figures 4 and 5 together display the results of this activity. In Figure 4, 
the fluid flow processes and their interactions are declared and in Figure 5 the 
thermal processes and their interactions are displayed. It is important to note that 
they axe not separate, they interact significantly. It is only for convenience in the 
presentation that they are placed in different figures. 


5.4 - THE PRIMITIVE EQUATIONS. 

The modeling system provides a basis for ordering the equations and 
labeling the variables of the model. First, for the matter exchange and storage 
portion of the modeling system, and then for the thermal exchange and storage 
portion of the system, we will list the primitive equations in the sequential order of 
the point numbers and/or path number(s) to which they refer. Reference to Figures 
4 and 5 should be made as needed. The choice of representation of the flow 
processes via “compressible” or “incompressible” flow models is made based on the 
process being described. In the present case, the pressure variations across the 
modeled processes are rather small, and thus incompressible representations are 
appropriate. 

The procedure is the same for each of the processes. First, the process is 
identified and then the modeling process equation which is needed is stated, first in 
its typically standard form. In many cases, these expressions are nonlinear. Next, 
supplementary equations are presented in order to demonstrate how relations 
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among the variables may be brought in to convert the expressions to relate only 
the selected set of variables. After this is done, the individual expressions are 
linearized about the prescribed operating condition. Then we move on to the next 
process or conservation constraint. In some cases, a process which might very well 
be included is stated, but because of the early phase nature of this study, 
insufficient physical configuration definition has been done, and it is judged that a 
meaningful study can be conducted though it is not included. “Meaningful” is a 
quality that is in the mind of the investigator; here we mean that we judge that a 
majority of the knowledgeable persons in the field would agree that it was a 
reasonable and understandable thing to do and that a worthwhile study could be 
conducted without accounting for the process. 

5.4.1 - SOME BASIC EQUATIONS. 

Some frequently used models of commonly encountered physical processes 
are listed for easy reference. 

Ideal gas law: PjVj = Nj Ft T- 

First law of thermodynamics for steady state, steady flow: h out - h in = Q + W 

N i 

Molar density: p- = rA 

J j 

Molar flow rate: n, = V i p, A, 

Total pressure loss model: ideal gas, steady, adiabatic, no- work flow: flow from 
point 3 to point 1 

( P 3 + \Pm z V) - ( p l+ Ki V l 2) = AP 3-i 
4 

n t = n, k ; Jt = 1 => nitrogen, fc = 2 => oxygen, k = 3 => carbon dioxide, and k = 4=> 

k = 1 

water vapor 

Molar content of the j.— volume: 
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N j 


4 

£n. k 

J’ 


fc= 1 


Molar concentration: of the 


</» 


N 


Cc j,k= 


j.* 


N i 


constituent in the f h volume: 


5.4.2 - THE MATHEMATICAL MODEL OF THE INDIVIDUAL PHYSICAL 
PROCESSES IN THE MODELING SYSTEM AND OF THE CONSERVATION 
CONSTRAINTS 


Fluid inert ance of the air column in the upper duct system (path i): 


1000 gc A t 

n i “ MWjL^ * P 35 p 3l 

^1 = ^23 

with linearization: 

_ 1000 gc £D 

6n i - MW.tr ( 6P 35 " * P 1> 


Rate of accumulation of the air constituents in the shoot upper chamber (point l)i 


Using the relations 


N 1 it 

n 3 , k ~ n 3’ k ~ X ’ 2,3 ’ 4 


ni ’ fc_ N 36 

n 14 = n l + n 18 

one obtains: 
Nitrogen. 


N 36,Jt . _ N 36,k _ _ , _ , 

n 14 t k — 1,2,3 ; ~ n i4 n 18’ k ~ 4 


36 


N l,i ~ n i,i ~ n 3, l 
• ^36 l ^1 i 

N 1,1 = ("1 + n is) - n 3 ITT 


with linearization: 


18 


N 6No fi Nog 4 

= (<n, + <n 18 ) + K + n 18 ) - K + n 18 ) ^^36.* 


— 6 no 


N, , «N, . No 4 

^i-n 3 ^ + " 3 ^E/ N M 


Oxygen. 


N u = n 12 - n 2>2 [H - n 2 2 [resp] - n 3 2 


N, 


No 


i og 2 1 1 1 2 

N lf2 = (nj + n 18 ) - n 2 2 [H - n 2>2 [resp] - n 3 

36 ^ 

with linearization: 


N 36,2 ± 


«*!., = (in, + ft.,,) ^ + (», + »,s) ^ - (». + ”.8> ^ £ <N 3M 


No 7 N 1 2 4 

Caxbon dioxide. 

^1,3 = n 1.3 ~ n 2,3(H “ n 2.3t res P] ~ n 3,3 


N 


36 3 j, 

^1,3 = ( n l n 18) ^Sj -1 n 2,3l^^l ” n 2,3[ reS P] ” n 3 “N“ 

36 1 


N 


1,3 


with linearization: 


6 N 


N 36,3 , , iN 36,3 


N 36,3 4 


"l,3 = ^"1 + 6n is) T^J- + ( n i + n is) ~W^r ~ (n i + n i8 } l^ fc 5/ N 36,fc 


36 


'36 


— 6 n 


N 1,3 _ «*1,3 . _ N l,3 < 


3 iff - "3 -Rf + "3 jf 


Water. 


^1,4 _ n l,4 - n 2, 4 n 3,4 


^ . , N 36,o . N l,o 

N 1.0 - ( n l + n 18) " n 18 * n !,4 n 3 loj - 

with linearization: 
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N 




36,4 




1,4 — ^ n l in is) N ' Dl8 ^ 

36 


36,4 


N 


36 


- K + "i 8 ) 


N 


36,4 


n 36 2 *= l 


£* N 36,fc “ ^ n i8 


N i « <N,. N, . 4 

- ^ n 3 ~ n3 “^7“ + n3 N^fcli^ 1 .* - 6ri2 ' 4 

Total molar content of the upper chamber air. 

"i = iX.- 

t = 1 

with linearization: 

<N l = i<N,, 

l = l 

Energy accumulation in the shoot side upper chamber. 


Ej = E A,.* - E h 3tk - E h, k + (UA) 2 ( t 2 - Tj) + (UA) 4 ( t 4 - t l ) - Pj v : 

1 k=l k=l h=2 

where 

(UA) 2 = h A pf 


h C2 = 11-30 ( 


*3 \°- S 


c 2 V 0.1p 2 (A p/ +A 3 ) 


(UA) 4 = h C4 A 4 


2 Nu t k. 


1. = 


Nu , = 0.332 Re ,/, Pr, 1/3 

*= L 1 Lj 1/2 1 


Pr 


1 ~ 



Re L 


1 


p 2 MWj V 1 Lj 
1000 JT { 
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Inserting equivalent thermodynamic properties; 


( N u + N lf2 c Vi 2 + Nj 3 c Vj 3 + N m c Vj 4 ) fj = 

( ° 1 ' 1 C P35,i + Bl ' 2 Cp 35,2 + ni - 3 Cp 35,3 + ni - 4 Cp 35,4 } (T 35 " 273 ’ 15) 

- < n 3,l c Pl J + "3,2 c Pl 2 + n 3,3 c Pl 3 + "3,4 c P j J (Tj - 273.15) 

" " 2 . 2 PH c p 2 2 (T 2 - 273.15) - n 2 2 [resp] c pj ^ (Tj - 273.15) 

- »2,3[H c Pl 3 (Tj - 273.15) - n 2 3 [resp] c p2 * (T 2 - 273.15) 

- "2,4 ‘p 2 4 (T 2 - 273.15) + h C2 2 LAI A p/ ( T 2 - Tj) + h^A, ( T 4 - Tj)- Pj Vj 

Using 

" lCp 35 ” Si DhkCp 35,k 

" 3Cp l " fc ?i n3 ’ tCp l,it 
and linearizing gives: 

* N U Cv l,i + N h2 c vj 2 + n 1, 3 c v l 3 + N i, 4 c v lt4 ) = 

n ! c P 35^ T 35 + ^ n i c p 35 ( T 35 ~ 273 - 15 ) ~ "3 c Pl ^ T l + ^"3 c Pj( T l ~ 273.15) 

- "2, 2 [H Cp 2 2 6T 2 - Sn 2 Jhv] c p2 2 (T 2 - 273.15) 

- n 2 2 [resp] c pj ^ 6Tj - 6n 22 [resp] c pj ^ (Tj - 273.15) 

- n 2 , 3 [H c Pl 3 - 6n 2 3 [H c pi g (Tj - 273.15) 

- " 2 , 3 [ res P] c p 2 3 *T 2 - «5n 2 3 [resp] c p2 3 (T 2 - 273.15) 

- "2,4 c p 2 4 6T 2 - 6n 2 4 c p2 ^ (T 2 - 273.15) 

+ h C2 2 LAI A p/ ( 6T 2 - 6Tj) + h c< A 4 ( 6T 4 - 6T j)- Pj 6Vj 
Equation of state of the air in the shoot upper chamber. 

P, Vj = Nj R T, 
with linearization 


£ p l v i + p i^ v i = ^ N i RT i + N i R ^ T i 

Airflow total pressure loss: from shoot upper chamber to shoot lower chamber 
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with linearization: 


«*1 - 6P b = 


2 n 3 MW 3 
2000 Pl A 3 2 g c 1,1,3 


in-, 


N 


l.fc 


MW 3 = EMW,^ 


A 2 n 3 MW 3 

3 3 2000 pj A 3 2 g c 

5PJ _ 6P 5 = K 3 6n 3 


A 5 - Aj 


Rate of accumulation of the air constituents in the shoot lower chamber {j>oint jQi 
Using the relations 


N 


5, it 


n 5,*--N^r “s 

and 


n 5 , k = 1 , 2 , 3 , 4 


n 5 — n 9 "I" n 10, 2 

one can obtain the following expressions for the accumulation of material in the 
lower shoot growth chamber. 

Nitrogen. 

^5,1 = n 3.1 " n 5.1 


N 


l.i 


N 


5,1 

*5,1 = “3 TTj “5 ”N^~ 


N 1 l N 5 .i 

^5,1 = n 3 -JT” “ ( n 9 + n io, 2 ) T^~ 
with linearization: 


6 N 5 j — 6 n 3 


^ + ”3 ^ - (<n * + ^ k 


iN 5 1 N 5 ,i 

- (n 9 + n 10 2 ) -jjt“ + ( n 9 + " 10 , 2 ) TfT .?/ N 5 ,k 
5 5 
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Oxygen. 


N -n ^ ^ 

*5,2 - n 3 Nj n 5 N 5 

*5,2 ~ n 3, 2 - n 5, 2 

N 1>2 N 52 

N 5,2 — n 3 “n” ( n 9 + n 10, 2 ) “NT - 

1 5 

with linearization: 


N, 01 N, n, 9 4 

**5.2 = *"3 -NT- + n 3 -NT 2 “ 113 FT E *Nj f * - (fa. + fa 10t2 ) 
1 1 2 fc=i 


<N 1.2 . N l ,2 ± 


2 Nt 2 4 

- (n 9 + n 10 2 ) h (n 9 + n 10?2 ) 


Carbon dioxide: 


N 


1,3 n 5,3 


N c 


'5,3 - "3 “N^ n 5 lsr 5 

N 5 ,3 = n 3,3 “ n 5,3 

N - n Nl ’ 3 ( \ *5,3 

*5,3 ~ n 3 V ( n 9 + n 10, 2) "M~ 

1 5 

with linearization: 


**5,3 = 6n 3 


N 


1,3 


S N 


1,3 


*1.3 A 


Nf +n 3 1^- n 3 ^ E™, * - (fa. + fa 10, 2 ) 


6 N 


( n 9 + n 10,2) N 5 5 3 + ( n 9 + n 10,2) •j^T j £^ N 5 <k 

Water. 


N 5,3 A 


*5,4 ~ n 3,4 - n 5,4 


N 1 . Nr 

j 1.4 5,4 

*5,4 n 3 Nj 


N 


1,4 


N 


*5,4 ~ n 3 ( n 9 + n io, 2) “N 

with linearization: 


5^4 

5 


**5,4 = Sn 3 


N 


1,4 


**1,4 _ *1,4 < 


*1 + " 3 ""nT" - 113 NTs N l- fc _ ^" 9 + * ni0 - 2 ) 


<5N 


( n 9 + n 10, 2) W + ( n 9 + n 10, 2) v ’2 £**5 Jt 
n 5 fc=i °’ K 

Total molar content of the shoot side lower chamber: 


*5,4 < 



The total molar content of the shoot side lower chamber is given by 

"5= IX. 

i = 1 

with linearization: 

<N 5=X>5., 

1 = 1 

Energy accumulation in the shoot side lower chamber . 

e 5 = E - t i,k + (ua ) 7 ( t 7 - t 5 ) 

k=i fc=i 

1 

< N 5.i c v 5i1 + N 5.2 c v 5 2 + N 5.3 c v 5 ,3 + Cv 5,4 } T 5 = 

( n 3,l c Pl , + n 3, 2 c Pl 2 + n 3,3 c Pl 3 + n 3,4 C P j J ( T 1 “ 273 ’ 15 ) 

- ( n 5,l C P5 , + "5.2 C f ) 5 2 + n 5, 3 C P5 3 + n 5,4 C p 5 J ( T 5 " 273 15) 

+ ua 7 ( t 7 - t 5 ) 

Including the simplification 

( n 9 + n 10,2) C P& - n 5,l C p 5 J + n 5,2 C p 5 2 + n 5,3 C J > 5 3 + n 5,4 C J >5 4 

along with linearization gives: 

( N 5.i Cv 5,i + Ns - 2 Cv 5.2 + N5 - 3 Cv 5,3 + N& ’ 4 Cv 5,4 ) = 

n 3 c Pl 5T l + ^ n 3 c Pl ( T l “ 273 ' 15 ) “ ( n 9 + n 10, 2 ) C P 5 5T 5 

- («n 9 + fa 10i2 ) c p5 (T 5 - 273.15) + UA 7 ( 6T-j - 6TJ 
MW 5 = MW! 

(UA) 7 = h c? a 7 



Nu T = 0.332 R*x 1/2 P t * 1 ' 3 
x=L 5 l 5 a 

p 5 mw 5 p 5 l 5 

Ke L 5 _ 1000 p 5 

L 5 = ^chamber ~ ( Kanopy + d PPFre/) “ d P/ 
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Equation of state of the air in the shoot lower chamber 

P 5 V 5 = N 5 RT 5 

with linearization: 

* P 5 V 5 + I V V 5 = * N 5 RT 5 + N 5 R * T 5 

Airflow total pressure loss: flow from lower chamber into discharge duct (path 5). 
MW fi RTr , 

P - P —Y 6 5 1 2 

r 5 r 6 _1V 5 2000 P 5 a 5 2 6 

with linearization: 

MW 6 RTc- , 

- 6Pg =K 5 — p 2 - — L n 6 <5n 6 (the incompressible flow model of this process is 

r 5 A 5 

employed) 

MW g = MW 5 


Pe — Ph 


Airflow total pressure loss: flow duct (path 5 between points 6 and 8); both 
direction change (“bend”) and duct length. Both Ignored: P g = P g 

Airflow total pressure loss: in expansion between points 8 and 9 (path 5 ) 

Ignored: (P g = P g ) 

Airflow total pressure loss: HEPA filter (path 6) 

P 9 “ p 10 = K e n 5 

K 6 is essentially a porous plug characteristic. Linearization gives 
6P g — ^P^q =2A" 6 n 5 6n 5 

Airflow total pressure loss: from HEPA filter discharge through control valve 12 
(path 7). Includes contraction loss in reentering duct upon exiting HEPA filter, the 
fixed (wide open) loss of the valve, and the loss due to the valve element being 
partially closed. A„ is the control variable. 

P 11 - P 12 = n s 2 
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k 7 = k$+k\ 


k° 7 = 


K? MW 7 RT n 


4= 


kJ mw 7 rt h 



linearization gives: 

— ^12 ^7 n 5 2 + 2^7 ^5^5 

MW 7 = mw 5 

All the pressure loss from point 5 to point 12 can be combined. 

P 5 - P 12 = i K s + K 6 + ^7 + ^ 7 ) n 5 2 

linearization gives: 

$P 5 - «P 12 = 2 (/f 5 +tf 6 + tfrXng + OlO,2)( 5n 9 + 6n 10,2) + ( n 9 + n 10,2) 2 6K 7 


Oxyg en extraction: at point 12. 
n 8 = n 5 — n io, 2 

Airflow total pressure loss: due to extraction of oxygen at point 12. 


Ignored. 

Airflow total pressure loss due to flow shear at duct wall: all the lower duct system 


friction-induced pressure loss is lumped into one single section in this model (path 
13 between points 12 and 13). 
p 12 “ P 13 = n * 2 

* - <.. 7 3 5 )(1^)(IV\^> 

n 5 = n 6 = n 7 = n 8 + n 10,2 

Ps MWg P 8 Lg 


Re. — 
"8 


1000 /i 8 
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V = _^8_ 

8 Ps ^8 
Ps ~ P% 

MW 8 = EMW U -A* - mw 7 , 2 ^ 

n 9 = n 8 

MW 8 = MW 9 
linearization gives: 


^ P 12 - ^ P 13 = 2 K 8 n 9 in 9 


Fluid inertance in the lower duct: 
_ 1000 g c A 9 


" 9 ~ MW 9 L 9 

^9 = ^8 


Lq ~ L 8 


( p !3- p i4) 


14' 


MW 9 = MWg 
linearization gives: 


Air Mover and Drive Motor Subsystem (AMfcDMSl. 

It was concluded that it was desirable for the model of the AM&DMS to 
reflect the addition of power to the fluid-thermodynamic systems and the rotational 
dynamics of the rotor components. To accomplish this, the relations among 
pressure rise, torque, rotational speed, molar flow rate, power transfer, and the 
rotational dynamics equations must be generated. 


The Fan Equations. 

An option to using First Principles in the modeling of the air mover 
characteristics is the use of the empirical fan equations which are implicit in the 
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conventional fan curves so common in industry; 


P 16 _ P 15 - f /an ( w 37’ n n) 
linearization gives: 


12 


m f „ m jan ( W 37' n n) t 9f /on ( w 37- n 12) , 

6p ig ~ 6p ib = K ^ 7 6w37 + dKT 2 6n 

this is a steady flow relation, and approximations of the partial derivatives can be 
taken directly from fan curves. 

Power Conversion (mechanical rotational power to fluid power): 


With perfect power conversion, the mechanical rotational power supplied to 
the fluid via the air mover impeller would be converted to fluid power. Since the 
process is not perfect, the overall efficiency is used to describe the actual situation. 
Mechanical power in (delivered to blower rotor) x efficiency = increase in fluid 
power out 

n 1 2 

u 37 T 36 % ~ ( P 16 ~ P 15^ Pj 4 

assuming constant density, linearization gives: 

w 37 Sr 36 ? o + 6w 37 r 36 % = ( P 16 " P 15^ ~T + < iP 16 " iP 15) ~T 

Air Temperature Rise Due to Air Mover. 


The temperature rise associated with the passage of the air through the air 
mover is the sum of that due to the isentropic pressure rise: 

n 1 2 

n i2 c pi2 t 16 - n i2 c pi2 t 15 = ( P 16 _ P 15^ 

and that due to the overall inefficiency of the air mover: 

n 12 c pl2 T 16 " n 12 c pl2 T 15 = ^ “ 1o) u 37 t 36 
Linearization of each (assuming constant density) gives: 

c pi2 6t 16 ~ c p ^ T 15 = ( ^ P 16 “ ^ P 15^ pJl 
and 

c pl2 ^ T 16 “ c p 15 = 0 - 37 r 36 + (! - 37 ^ r 36 
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n 12 ~ n 9 

Combined pump and motor rotor dynamics. 

Drive motor rotor and pump rotor: 

(J impeller + ^motor) ^37 + ^l w 37 + ^2 W 37 = T c ~ r 36 

linearization gives: 

(J»mpe//er + ^motor) ^37 + (^ 1+2 ^2 W 37)^ W 37 = ^ T c ~ ^ r 36 

The Gas Injection and Mixing Chamber 

The gas injection and mixing chamber is modeled to serve two purposes: 
account for the injection of carbon dioxide, oxygen, water vapor, and nitrogen 
content into the flowstream, and to reflect the effects of the volume of the ducting 
and fittings in that area on the overall system dynamics. 


Airflow total pressure loss: into gas injection and mixing chamber. 


P 16 “ P 36 “ 


k 13 mw 13 rt 16 


2 P 16 8c A 


13 


0-^)n 13 

a 14 


mw I3 = mw 9 

this pressure loss is ignored 

Rate of accumulation of the air constituents in the gas injection and mixing 
chamber (point 36): 

In the following, the equations of like level are grouped. Previously all 
equations addressing the same constituent were grouped. 

^36,1 = n 13, 1 + n 17,l“ n 14, 1 
^36,2 = n 13, 2 + n X7, 2 — n 14, 2 
^36,3 = n 13,3 + n 17,3 — n 14, 3 
^36,4 = n 13, 4 + n 17,4~ n 14,4 
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Nr , Ng k 

n 13,fc = ( n 9 + n 10,2)’ k = 1,3 ’ 4; n 13 , k = "NT - ^ n 9 + n 10,2) _ n 10,2’ k ~ 2 


N 36,fc , v N 36,fc , , „ 

I,14,it = ~ PT 36 ~ " 14 - iq 36 ~ 1 + n is )> k ~ 1 ' 2 ' 3 ’ 4 


Nr 


N, 


, " 5, 1 , , V ‘'36,1 

^36,1 ” ( n 9 n 10, 2) ^ + n 17,l v n l“*" n 18) Ngg 


N ^36 

^36,2 = ( n 9 + n 10i2^ - n 10,2+ n 17,2 “ ( n l+ n is) -fl^T 


x N 5,3 , ^ v N 36,3 

^36,3 — ( n 9 + n 10,2) Ng n 17,3 ' n l + n 18) Ngg 


Nr 


Nr 


x '5,4 , , x 36,4 

^ 36,4 — ^ n 9 n 10 , 2 ) Ng + ni7 ' 4 n 18 ^ Ngg 

linearization gives: 

^ 36,1 = ^ n 9 + ^ n io> 2 ) ~^T + ( n 9 + " 10 , 2 ) U 7 ( n 9 + " 10 , 2 ) N 2 fc 5/ N 5 .*= 


N 5,i ^ 


n 36, 1 , N 36.- ^ 


+ ^ n i7,i _ (^ n i+ ^ n is) Ngg ^" 1+ ni8 ^ ,Nifi + n * 8 ^ Nflfi 2 ^^36,* 


N 


6N 


36 


^36,2 — ^ n 9 ^ n 10,2) Nc ^ n ® ni0 ’ 2 ^ Nr. ^ n9 ni0 ' 2 ^ Nr 2 


X, 


- £n 10 2 + £n 17 2 - (^nj+ 6 n 18 ) N ^_ ( n i + n 18 ) + ( n i+ n is) N ^2 ^ SN 36 ,* 


N 36,2 4 


'36 


Nc 6Ng 3 Ng 3 4 

^36,3 = + in 10,2) T^~ + ( n 9 + n 10,2) “N^ ("9 + n l<u) N ^2 

+ *n 17>3 - (in,+ Sn u ) ^ - (n,+ n 18 ) + (n,+ n 18 ) 

Nc 6Ng Ng 4 4 

^36,4 = ^ n 9 + ^ n 10,2) + ( n 9 + " 10 , 2 ) “ITT ^" 9 + ni0 ' 2 ^ Nr' 2 


+ ^ n 17,4“ (^ n l + ^ n is) Ngg ( n l+ n 18) Ngg + ( n l+ n is) N 3g 2 fc ?j^ N 3M 
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Energy balance in the gas injection and mixing chamber 


E 36 - .2 A 13 ,k + J2 h°7,k ~ A 14,k 

1 


fc=i 


k=i 


Jfc=l 


( N 36.i < ' v 36 j1 + N 36,2 c v 36 2 + N 36,3 c v 36 3 + N 36,4 c v 36>4 ) T 36 “ 

( " 9 - 1 c Pi6,i + n9 ’ 2 c pi 6,2 + n9 - 3 c pi 6, 3 + n9 - 4 (Tl « " 273 15) 

+ ( n 17,l C Pi 71 + "l7, 2 C P 17 2 + ni7 - 3 C Pl7,3 + "l 7 * 4 C Pl7,4^ ^ Tl7 “ 273 ’ 15 ) 

- ( "l4,l C p 36 j + "14,2 C p 3g 2 + "14,3 C p 3 g 3 + "14,4 C p 36 J ( T 36 “ 273 ‘ 15 ) 

n 9 C Pi6 “ n 9’ 1 C P16,1 + n9 ' 2 C Pl6,2 + " 9 ’ 3 C Pl6,3 + n9 ’ 4 C Pl6.4 

n l 4 c p 3 g = "14,1 C p 36 , + "14,2 C p 36 2 + "14,3 C p 3g 3 + "14,4 C p 36 >4 

linearization gives: 

( Cv 36,i + N 36,2 c v 36 2 + N 36,3 c v 3g 3 + N 36,4 c v 36 J * T 36 = 

^"9 c p 16 ( T 16 “ 273 - 15 ) + "9 c p 16 

+ ( *"17,1 c p 17 J + *>17,2 C Pj7 2 + 5n 17.3 C Pl7,3 + ^l 7 - 4 ^17,4) (T 1 7 ~ 273 ' 15 ) 

+ ( "l7, 1 C p 171 + n 17,2 C Pi7 2 + n l 7 t 3 C Pl7,3 + "l 7 - 4 C Pl7,4^ ^ 17 

— (^n a + ^n 18 ) c p 3 g(^36 - 273.15) — (nj+ n 18 ) Cp^6T 3 g 
Equation of state for the gas injection and mixing chamber. 

P 36 V 36 = N 36 R T 36 

iP 36 V 36 = ^36,*) R T 36 + N 36 R iT 36 

k = 1 

Airflow total pressure loss: exiting from the gas injection and mixing chamber. 

(GIMQ - 

MW, fi RTo fi , A, 2 

P _ P —V 36 3o 1 / a 14 ti 1 \ 2 

*36 M7 -K i4 2000 Pqa A 2 a 2 1 14 

36 a 14u a 14 d 

where 


A 14u is the flow area at the upstream end of flow path 14 
A 14J is the flow area at the downstream end of flow path 14 
MW 36 =£p 3 g ■ c p •; the molecular weight of the gas in the GIMC 
this pressure loss is ignored 
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Water Vapor Removal System 


The water removal system must route the necessary fraction of the total 
flow through the condensing heat exchanger so that the evapo-transpiration water 
is removed. It is presumed that the temperature of the air exiting the condensing 
heat exchanger is fixed (e.g., by the use of feedback control) at a prescribed level, 
and is saturated at that temperature. 

By-Pass Fraction Determination 

The by-pass fraction is the fraction of the total flow through the water vapor 
removal system which is diverted around the heat exchanger 

A n 16 _ n 14 - n 15 
Xi 4 - n 14 - n 14 

There are three equivalent, relevant, and related relationships among specific 

humidity, concentration, and partial pressure, 

(1) SH 14 = Ngg J ^ Njg (2) SH 14 = Cc 14 4 /(1 — Cc 14 4 ), and 
1 = 1 

(3) Cc 14 4 = (partial pressure of the water vapor at 36)/Pjg 

The water removed by the heat exchanger can be expressed in terms of the total 

air side flow by 

n i 8 = (SH 36 - SH 21 XI - SHjg) n 15 

and linearization gives: 

6 n 18 = (^SH 36 - tfSH 21 )(l - SH 36 ) n 15 - (SH 3g - SH 21 ) 6 SH 3 g n 15 
+ (SH 3 g - SH 21 )(1 - SH 3 g) 6 n 15 

n 15 = (1 - X 14 ) n l4 

linecirization gives 
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6 n l 5 = - 6 Xu n 14 + 0 - X 14 ) <S 4 

The nominal value of Xj 4 is established by setting n 18 to equal n 24 to determine the 
nominal value of nj 5 , and then solving for Xi 4 - The requirement of matter 
conservation in the flow split establishes that 

n 15 + n 16 = n 14 = n l + n 18 

which with linearization gives: 
feu + in 16 = in! + in 18 

Airflow total pressure loss: through control valve 21 (heat exchanger by-pass) 

P 18 = P 17 
P 19 = P 17 

as is pressure loss in the by-pass up to the control valve. 

P 20 = P 18 


P 20 ~ P 23 “ 


K2qMW 36 RT 36 


2 P 


O-^f 2 ) n : 


20 g c A 


t/20 


A 2 

a 20 


20 


P 20 “ P 23 ~ *20 n 20 


K 20 MW, fi RTo 8 A 2 ™ 

IT -I- 20 36 00/1 ^v20\ 

*20 “777 7~T2 V i “-72“V 


2 Pon g c A 2 


20 Sc A t>20 

MW 20 = MW 36 
linearization gives: 


A 2 

a 20 


^ P 20 ^ P 23 — ^ n 20 *20 ^ n 20 + n 20 2 ( -’ A : 


20 


n 21 = n 20 

Airflow total pressure loss: through control valve 22 and dehumidifying heat 
exchanger. 

Control valve in water removal system heat exchanger path duct: 

P 21 _ P 22 = *19 n 2i 2 
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. k 19^ W 36^ T 36/i A «19'| 

*19 =~^ — TT2 ' A 2 ’ 

2 ”21 g c A «19 A 19 

linearization gives: 

6P 21 - SP 22 = 8K l9 n 21 2 + 2tf 19 n 2) 6n 21 
For the heat exchanger: 

P 19 _ p 2 i = ni9 

For preliminary work, one can estimate K 1S via Reynolds’ analogy or simply budget 
a reasonable pressure loss. Linearization gives. 


6P 19 - * P 21 = 2A 1S n 19 6n 19 

These two pressure loss expressions can be combined: 
*P 19 - 6P99 = n 19 2 6tf 19 + 2(^5 + K l9 ) n 19 6n 19 


22 “ “19 ^19 


mw 22 = mw 21 = mw 14 

-4- N 36,fc 


mw 22 =Emw 36i11 


- MW 4 JJ7T 


Jb=l 

n 22 = n 15 — n 18 


36 


15 


Energy balance: condensing heat exchanger, (reference is liquid water at triple 
point) 

“is c p, 36 T 36 “ “22 c p , 21 T 21 ~ n i8 A 38 = 0 
linearization gives: 

“> 5 %,36 T 36 +". 5 'p, 36 <T 36 %,21% %,21* T 21 - <n » *38 + “>*“38 =° 

Energy balance: recombination of flow streams of the water vapor removal system. 

“l c p 24 ( T 24 - 273 15 ) = n 2i C P 2 3 (T 23 ' 2?315) + “ 22 Cp 22 (T22 ” 2?3 * 15) 

4 

“i C P24 = fc 5 1 " 1 ’ fc Cp 24, fc 
4 

n 2l Cp 23 - Cp 23,fc 


"22 Cp 22 - t 5 1 “22. fcCp 22.fc 


n 22 C I 

n 2 i = n 20 = n 2 - n 10 2 - n 15 


Cp 23 = Cp 36 
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Cp 22 = Cp 21 

Cp 24 =Cp 36 = Cp 21 
T 23 = T 36 

T 22 = T 21 

linearization gives: 

n i C P 24 iT 24 + 6n i c p 24 ( T 24 ~ 273 ’ 15 ) = ( n i " n io ,2 - n is) c p 36 ^T 36 

+ (6n, - 6d 10 2 - 6n 15 ) (T 36 - 273.15) + (£n, 5 - 6n 18 ) c p ^ (T 21 - 273.15) 

n i 6T 24 4- 6n, (T 24 - 273.15) = (n, - n 10<2 - n 15 ) 6T 3fi 
+ (<511! - <5n 10 2 - <5n 15 ) (T 36 - 273.15) + (<5n, s - 6n, 8 ) (T 2J - 273.15) 

Airflow total pressure loss: combining flows of water removal system paths 

(1) heat exchanger bypass path (path 16) 

P 24 - P 23 = ~~~~~ ^3 ^^ RT 25 j ( ^3 )2 _ (£1)2 _ ( £^) _ ( ^_ M 

24 A 23 A 21 A 23 A 21 A 22 A 23 

This pressure loss is ignored. 

(2) heat exchanger path (path 15) 

K 22 73 n n n 2 2 

P«,. - P„„ = — — — 22 21 r , /_23\2 „ / n 22\2 _ , n 22 \ / n 22 m 


P 24 " P 22 “ 




This pressure loss is ignored. 

Airflow total pressure loss: duct between water removal subsystem and air 
temperature manipulation subsystem (path 23). 

P 24 “ P 25 = ^23 n 23 2 

K a = (.0735)(^iiS)(R« L )4 ( MW 23 , 

A 23 l 23 7 v 2000 p 23 A 23 2 


d- _ P 23 MW 23 ^23 ^23 _ n 

L 23 " 1000 /i 23 ’ V 23 - 

MW 23 = MW, 


“23 

p 23 A 23 


linearization gives 

^ P 24 “ ^ P 25 = 2 K 23 n i ^ n i 
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The Airflow Cooling/Heating System 

The airflow cooling/heating system (AC/HS) must route the necessary 
fraction of the total flow through either the cooling heat exchanger or the heating 
heat exchanger so that the temperature of the flowstream exiting the AC/HS is at 
the required level. It is presumed that the temperature of the air exiting the active 
heat exchanger is fixed (e.g., by the use of feedback control) at a prescribed level, 

and no condensation occurs. 

By-Pass Fraction Determination 

For the present demonstration, the by-pass fraction is the fraction of the 

total flow which is diverted around the cooling heat exchanger 

A n 25 _ n 23 ~ n 26 
X23 1*23 **23 

n 25 = *23 n 23 

linearization gives: 

^ n 25 = ^*23 n 23 + *23 ^ n 23 

The nominal value of x 23 is established by performing a thermal energy balance 
analysis so that the required nominal discharge temperature at 35 is achieved, as is 
explained in the following. 

Matter conservation across the flow split requires that 
n 25 + n 26 = n 23 

and linearization (in this case, a redundant operation) gives: 

^ n 25 + ^ n 26 = ^ n 23 
n l = n 23 

Energy balance at recombination of airflow st£eams from airflow cooling systegk 

n i c p 35 ( T 35 " 273 15 ) = n 3i C P 33 < T 33 " 273 15) + n 32 C P34 (T34 " 273 ' 15) 

+ "30 C P30 (T32 - 273 - 15 ) 
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where 


" lCp 35 “ £ 1 Rl ' k %5,k 

4 

n 31 c p 33 ~ fc E n 31,fc C P 33 jb 


1132 Cp 34 " fc 5i n32 ’ tCp 34,* 

n 31 “ n 28 
n 32 = n 27 

Cp 33 = Cp 24 

Cp 34 = Cp 29 
T 33 = T 24 

T 34 = T 29 

along with linearization, gives: 

n i Cp 35 6T 35 + 6n i C P 35 ( T 35 - 273 15 ) = n 28 ^p 24 ^ 24 + *n 28 c p ^ (T 24 - 273.15) 
+ (Si h - 6n 28 ) c p ^ (T 29 - 273.15) 
or, since 

Cp 35 “ Cp 24 “ Cp 29 

n i * T 35 + ^ n i ( T 35 - 273.15) = n 28 6T 24 + 6n 28 (T 24 - 273.15) 

+ (6nj - in 28 ) (T 29 - 273.15) 

MW 26 = MW 27 = MWj 

P 28 = P 25 

P 34 = P 35 


One may establish the nominal value of X23 f° r the scenario demonstrated in 
this study through simultaneous solution of the two equations 
n 23 c p 35 ( T 35 “ 273.15) = n 25 Cp^ (T 33 - 273.15) + n 25 c p ^ (T 34 - 273.15) 
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and 


n 25 — *23 n 23 

to satisfy the requirement for T35 

Airflow total pressure loss: air entering each of the three airflow temperature 
manipulation subsystem ducts; 

(1) heating heat exchanger path (path 24). 

K 2 3_24 MW 23 RT25 


P 25 “ P 26 “ 


25 


[ (^) 2 - a, (^) 2 - a, (^-)) 

a 24 "23 


A 2 4 A 23 


This pressure loss is ignored 
(2) bypass path (path 25) 


P 25 “ P 27 “ 


^23 - 25 MW23RT25 


25 




x 23 


25 a 23 


)] 


This pressure loss is ignored. 

(3) cooling heat exchanger path (path 26) 

„ „ *23-26 MWl 3 RT 25 _ n 26 2 ^ n 23 x2 _ B ( n 26 p 23 M 

P 25 " P 28 = l a i <A^ “ ( A^> ^ ( A 26 A 23 > J 


25 


*23 


Airflow total pressure loss: air flowing in each of the three airflow temperature 

manipulation subsystem ducts; 

(1) heating heat exchanger path (path 24). 

K 24 MW?, RT<; 


P 26 “ P 31 “ 


' 23 
3 26 


45(^ )2 

a 24 


This pressure loss is ignored. 
(2) bypass path (path 25) 

_ 4 f 25 L 25 MW 25 RT 25 


P 27 “ P 30 ~ 


D 


25 


2 P 27 A 25 


2 n 25 


This pressure loss is ignored. 

(3) cooling heat exchanger path (path 26). 

K 26 MW 23 RT25 ,n 26x2 

P VT — ) 

*28 


P 28 ~ P 29 “ 
linearization gives 


v 26 


iPjg - ^29 " ^26 ^26 

this pressure loss is combined with the fixed part of the control valve pressure loss. 
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Airflow total pressure loss: control valves in the air temperature manipulation 
system 

(1) heating heat exchanger path and control valve (path 29). 


heat exchanger; use Reynolds’ analogy or budget a reasonable pressure loss 

_ MW„ RT 0 
P 26 -P 31“ OD . 

2 P 26 g c A 24 

linearization gives: 


2 n 29~ ~ ^29 n 29 


^26 “^31 “ ^^29 n 29 ( ^ I1 29 


control valve: 


P 31 “ P 32 “ 


k 29 mw 23 rt 31 


2 P3I 8 C A „29 

linearization gives: 


A 

A 


(1 — *2 9 ) n 29 2 “ ^29 n 29 2 


29 


^ P 31 ” ^ P 32 “ ^^29 n 29 2 + 2A' 29 n 29 Sn 


29 


This pressure loss was not relevant to the scenario demonstrated in this study. 

(2) bypass path (path 28). 


control valve: 


P 30 “ P 33 “ 


K 28 MW„ RT 


2 p 3n 8c A 


“"'31(1-^5)112, 


30 ®c A »28 


v 28 


P 30 ~ P 33 - K 28 n 28 2 
linearization gives 


tfPjO ~ ^f*33 — 6^28 n 28 2 + 2^28 n 28 ^ n 28 

(3) cooling heat exchanger and control valve (path 27). 

heat exchanger: use Reynolds’ analogy or budget a reasonable pressure loss, 

P 28 — P 29 = ^26 n 26 2 
linearization gives 


6 p 2g — 6P29 — 2A" 26 n 26 6n 26 


control valve: 


P 29 “ P 34 “ 


K 27 mw 23 rt 31 


2 P 29 g c A 


(l-^)n 


27 


t>27 


v 27 
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P 29 -P 34“*27 n 27 2 

linearization gives: 

6 P 2 g — 6 P 34 = ^27 n 27 2 + 2 K 2 7 n 27 6n 2 7 


n 27 — n 26 


Airflow total pressure loss: recombination of air temperature manipulation 

subsystem flow (into path 1). 

(1) heating heat exchanger path (path 30) 


P P - K 25 - 33 MW l RT 35 r (^28 )2 _ - (^)2 _ ( "» . _ ) _ & ( x~^ ~)] 

Pqo - Pqr - P„ r 1 1 ^A 28 ^ 2 a,/ 3 ^A 23 A 21 ; A 22 A 23 


35 


32 “ r 35 
This pressure loss is ignored. 

(2) bypass path (path 31) 

^22-23 MW 22 RT2^ 


P 33 " P 35 “ 


22 




•)1 


This pressure loss is ignored. 

(3) cooling heat exchanger path (path 32) 

K 23 - 25 MW 23 RT 25 f „ ^ n 25\2 a ^^2 _ ^ ( n 25 n 23_^ 
P 34 ~ P 35 = A^j &2 ( A 23 } 43 A 25 A 23 ' J 

This pressure loss is ignored. 


MW 28 = MWj 

n l = n 27 n 28 

P 30 = P 25 
P 33 = P 35 

The heating heat exchanger control valve is fully closed in the scenario for which 
this study is done. The overall pressure loss for the airflow heating/ cooling 
subsystem can be expressed in terms of either of the two active paths ( 25 - 28-31 or 26 - 
27 - 32 ): 

6P 25 “ * P 35 = 2 ( n l “ n 2?) *28 ( 6n l - 6n 27) + ( n l - n 27) 2 6K 28 
SP 25 - ^P 35 = 2 (K 26 + K 27 ) n 27 6n 27 + 6K 27 n 27 2 
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Airflow total pressure loss from inside the discharge end of the flow duct (path 1) 

into the shoot upper chamber: 

MW, 2 RT. Al 

p 3 - p , = Kj - - 1 , 3 [ 1 - —3] n 2 

P 3 A 3 a 1 

linearization gives: 

„ £n « _ . MW, 2 RT, a\ MW , 2 RTo A\ 

<p 3 - sp i = 2 K ' 1 1 - tt i - K > shr 1 1 1 - if: 1 


This pressure loss is ignored. 


k 3 a 3 


'<5P, 


Recognizing the effect of glove port activity on upper chamber volume 
completes the inventory of relevant processes: 


V 1 = V , 


6Vj = V a 


5.5 - THE MODELED SCENARIO 

The control study is directed to a plausible research experimental condition. 
The defining parameters of this experiment are: 

(1) . Temperature of the air exiting the upper chamber, T t : 20 C. 

(2) . Relative humidity of the air exiting the upper chamber, RH X : 70%. 

(3) . Absolute pressure of the air in the upper chamber, 101,458 Pa. 

(4) . Temperature of the upper chamber wall, T 4 : 22 C. 

(5) . Source of PAR: on. 

(6) . Leaf area index: 30. 

(7) . Energy radiated from lamps which enters the CGC: 6000 Watts. 

(8) . Canopy absorbs all incident radiation. 

(9) . One-sixth of the absorbed energy evaporates transpiration water 

(10) . Five-sixths of the absorbed energy heats the air flowing through the CGC. 

The system performance and physical configuration requirements, the 
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parameters of the prescribed experiment, along with enough bridging assumptions, 
when substituted into the modeling equations, generate the nominal values of each 
of the flow variables and point variables of the CSSEDS model. 

Solving the above nonlinear equations for the equilibrium values of the 
system variables provides one with the values of all the variables which can be used 
in the linearized equations to determine the dynamical properties of the system for 
values of the system variables in the neighborhood of those of the experimental 

scenario. 

5.6 - NOMINAL VALUES OF THE POINT VARIABLES AND THE FLOW 
VARIABLES 

A p/ = 2.0 m 2 , A 1= = 0.1297 m 2 , A 3 = 0.64 m 2 , A 4 = 7.791 m 2 , A 5 = 0.1297 m 2 , 
A 7 = 6.2402 m 2 

bl = 0.0 Newton-meters/radian per second, 
b 2 = 0.0013 Newton- meters/ (radian per second) 2 , 

(c p )nom = 29.1531 Joules mol" 1 K~\ (c ?1 2 ) nom = 29.4908 Joules mol" 1 K _1 

( C pj’*) nom = 37.0476 Joules mol' 1 K~\ (c Pl ’ 4 ) nom = 33.7301 Joules mol" 1 K' 1 

The above values are also used for these properties at points 2, 17, and 36. 

(c Pl ) nom = 29.14 Joules mol' 1 K -1 

The above value is used for this property at points 5, 12, 16, 17, 21, 24, 35, and 36. 

(c v ) nom = 20.8387 Joules mol -1 K~\ (c Vl 2 )„om = 21.2055 Joules mol" 1 K' 1 

(cv 1,1 )" 0 = 28.7333 Joules mol' 1 K“\ (c v ’ ) nom = 25.4158 Joules mol' 1 K _1 

The above values are also used for these properties at points 5 and 36. 
d p/ = 0.12439 m 


Dg = 0.4064 m 
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— q 37? ^ — 40 .6 Pa/(radian per second), 

0U) 37 

^ w 37> ^g) = i 792 Pa/(mole per second) 
dn 9 

Khamber = °* 7576 m 


^canopy ^PPF rc y ^chamber 
( h c 2 )nom = 34 ‘ 577 Watts m " 2 K ~^ 

( h c,)nom = 4 ’ 239 Watts m " 2 K_1 
4 

( h c 7 )nom = 7 * 3196 Watts m " 2 K ‘ a 


^impeller = 4 * 468 k 8 

^ motor = 1328 k 8 m2 

( k l)nom = 0.0263 watts m -1 K -1 


K filter = .09555 n m -2 mole -2 sec 2 , 

K 3 = 1.0 n m mole sec , 

A 5 = 1.0 n m -2 mole -2 sec 2 , 

Kg — Kj ilier = .09555 n m -2 mole -2 sec 2 

K 7 ° — 0.00016 n m -2 mole -2 sec 2 , 
Kg = 0.002031 n m -2 mole -2 sec 2 

K 7 l = 0.055332 n m -2 mole -2 sec 2 

~ 68.4 n m -2 mole -2 sec 2 

K 21 — A a 20 = 158.8 n m -2 mole -2 sec 2 

K 22 = 114.0 n m -2 mole _2 sec 2 

K 23 — 158.8 n m -2 mole -2 sec 2 

K 26 = 0.0 n m -2 mole -2 sec 2 

K 27 = 0.1362 n m -2 mole -2 sec 2 

K 2 g — 1.215 n m -2 mole -2 sec 2 

K 2S = 1.215 n m -2 mole -2 sec 2 


LAI = 30 (unitless) 

Lj = 3.85 m L 5 = .2544 m L g = 3.85 m L 23 = 3.85 m 

(MW 1 ) nom = (MW 35 ) nom = 28.9700 g/mol, (MW 3 ) nom = 28.8277 g/mol, 

(MW g ) nom = 28.8277 gm mole" 1 (MW 9 )„ om = (MW 8 ) nom = 28.8277 gm mole" 1 

(MWjt,,, = 28.8277 gm mole' 1 (MW 5 ) nom = 28.8277 gm mole" 1 

K c p 35 )nom = ( n 3 c Pl )„om = (102.26)(29.67) = 3034.448 Watts K' 1 

( n i)nom = 102.2369 moles sec* 1 , ( n i,i )nom = 75.944 moles sec* 1 

( n i 2 )nom = 24.54225 moles sec -1 , ( n i t 3)nom = 0.122844 moles sec -1 

( n l,4)nom = L6278 moles sec_1 ' 

(n 2 2 [re* p])nom = 0.0 moles sec \ 

( n 2,2 )nom = ( n 2,2 ['“'Dnom = -0.00014699 moles sec -1 , 
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r—l 


(n 3 , 2 )n om = 24 ' 5420 moleS SeC_l - 
(n 3 , 4 )„ om = 1 ' 6505 m ° leS SeC_1 

(n 6 )„ om = ( n 3 )nom = 102 - 1592 m ° leS SCC ‘ 
(n 9 )nom = 102.1591 moles sec" 1 


= 0-0 moles sec' 1 , 

(n 2t3 [re.p)) nom = (n 2t3 ) nom = 0.00013542 moles sec' 1 , 

(n 2 . 4 )nom = -0.022236 moles sec" 1 
( n 3 )nom = 102.1592 moles sec' 1 
(n 3 c pi ) nora = (102.26)(29.67) = 3034.448 Watts K 
( n 3 ,i)nom = (ni.l)nom = 75.944 moles sec- 1 , 

(n 3t3 ) nom = 0.122712 moles sec' 1 , 

(n 5 )nom = (“aUn, = 102-1592 moles sec- 1 , 

(n 8 )„om = 102.1591 moles sec’ 1 
( n 10 , 2 )no»n = 0.00014699 moles sec' 1 

(n 14 )„om = (“i + “18 )non, = 102-2595 moles sec" 1 (n 15 ) nom = 2.0861 moles sec 

/ n \ - = (n 01 ) = 100.195636 moles sec -1 

( n i7 i)nom = 0-0 moles sec -1 ( n i7,2)nom = 0.0 moles sec 

(n, T *) = 0.00013542 moles sec -1 ( n l7,4)nom = ®’® moles sec 

rJ) = fni *) = 2.0861 moles sec” 1 

V n 19/nom v 1I 15^nom 

(“i 8 )nom = - (“2,4 )nom = 0 022236 moles sec- 1 
( n 2 i)nom = Mnom ~ (“ 22 )nom=100.195636 moles sec" 1 
( n 2 2 )nom = («l 5 ) Mm - (“w)nom= 2 - 063864 moles sec- 1 

(n 26 ) nom = (“ 27 ) nom = 76-64 moles sec- 1 (n 28 ) n<)m = 25.5969 moles sec’ 1 

(“2 2 )nom = (“isUm ' (“i 8 )nom= 2 -063864 moles sec- 1 


(N U tyj j + N U c Vi 2 + N U c Vi 3 + N lt4 Cv 1 4 ) nom = 1778.2525 Joules 


sec 


-1 


(NxU™ = 83.2528 moles, 


N 


iii) = .74266 unitless, 




N 


(-f^nom = 001200 UnitleSS ’ 

(N 5 ) nof » = 28-3055 moles, 

/ N 5,l ' _ / n 5.U _ Ah = 0.74266 unitless, 

(-fT-Jnom - nom \ n 3 >nom 

5 

f N 3,2 , _ / n 5.j \ _ /OliZ) = 0.239936 unitless, 

)nom V n 5 /nom 1 1^/nom 


(-jq-)nom = ' 24 U “ 5t,eSS ’ 
(^oom = 01614 UnltleSS ’ 
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N 


5,3 


i-tiflnam = (4f) 


"3 3 

v n 3 >nom 


N 


5,4 


( ^ = r JLi\ = 

V ^ fnom V n 5 /nom 


4 

( J ’ 4 ) 

V n 3 / nom 


( N 36)non» = 4L622 m ° leS ’ 


No/? . 

(~ )nom = -742663 unitless, 
36 

N 36 3 

)nom = -0012013 unitless, 
36 


Nu t — 35.4046 unitless, 
*= l 5 


= 0.0012 unitless, 
= 0.01614 unitless, 


f^ 36 ' 2 ) 

V N 36 n ° m 

Nqr „ 

( 36 ’ 4 ) 

V N 36 n ° m 


= .239987 unitless, 
= .0161401 unitless, 


(Pl)nom= 101,458 Pa, 

(P 3 )nom = 101,458 Pa 

(Pi - P 5 )no m = 8-84 Pa, 

(P 5 )no m = 101,449 Pa 

(P 5 - P 6 )no m = 214.765 Pa 

(P 6 )no ra = 101,234 Pa 

(P 9 )™n = 101,234 Pa, 

( P 9-Pl0)- ra = 999.2 Pa 

(Pio)ncm = 100,235 Pa, 

(Pll)nom = 100,235 Pa 

( p il - P 12)no ro = 599.5 Pa, 

(Pl 2 )nom = 99635 Pa 

( P 12 - P 13)nom = 212 Pa 

(P i3 )nom = 99,423 Pa 

( P 13 “ P 14)nom = 0 

(Pl 4 )»om = 99 ’ 423 Pa 

( P 14 - P 16)nom = -3850 Pa 

(P 36 )n 0m = 103,272 Pa 

( P 17)ncm = 103,272 Pa 

(P 24 )nom = 102,472 Pa 

(P 25 )no m = 102,260 Pa 

(P 3 6)nor„ = (Pl7)nor„= 103,272.1 Pa 

(P 35 )non, = 101,458 Pa 

( Pr l)nom = 0-707 unitless, 

( Pr 5)nom = 0 -707 unitless 

R=8.31441 Joules mole “ 'K ” 1 

(^L^nom = 42,600 unitless 

(Re^ ) nom = 4.413 x 10 6 unitless 
s 

( Re L c )nom = 14,329 unities 
5 

SH 21 = 0.008595, 

SHag = 0.016405 

Wiinam = 2 »3.15 K 

(T 2 )„ om = 295.15 K 

(T 3 ) nom = 291.4179 K 

(T 5 ) nom = 293.173 K 


45 


(TjyJnom — (^-^nom ~ (^^nom — 298.1 
(T 24 ) nom = 294.8 K 
(T 27 ) no „, = (T 28 ) nom = 294.8 K 
(T 30 )„ om = 294.9 K 
(T 36 ) nom = 298.1 K 

(UA) 4 nom = 33.03 watts K" 1 

( UA >7 nom = 45 - 676 WattS K_1 

Vj = 2.0 m 3 V 5 = 0.68 m 3 

(^l)nom = 0.00001983 kg m" 1 sec" 1 

* 7 p = -6 

(^)nom = 4L62 moleS m " 3 
(Xi 4 )„om = °- 0204 (unitless) 


K (T 21 ) nom = 278.2 K 

(T 25 )nom = (T 24 )„om = 2 9 4 -8 K 
(T 29 ) noro = 288.85 K 
(T 35 ) non , = 291.4179 K 


(UA) a nom = 4149.24 watts K 1 


V 36 = 2 Aj L 23 = .99869 m 3 
(n b ) nom = 0.00001983 kg m _1 sec 

(A*i )„ om = 41 -63 mol m" 3 

( w 37)nom = 189 rad 860-1 
(x 23 )nom = - 8134 (unitless) 
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5.7 - THE CONDENSED PRIMITIVE EQUATIONS 

To express the system of equations in the format of Equations (1) and (2), 
we assign r/,<r,u, d elements to the perturbed variables, 6 ( • ) as is shown in Table 1; 

As an aid to the reader, the inverse assignments are listed in Table 2. 

The result of using these assignments and inserting the nominal numerical 
values into the linearized equations is: 

»h = 1.279 <r 1 - 1.279 <r 2 ( 5 - 7 - 1 ) 

j) 2 = 0.7427^! - 0.3163 t/ 2 + 0.9127t? 3 + 0.9127tj 4 + 0.9127i/ 5 + 0.6330tj 13 - 1.826tj 14 - 1.826r? 15 - 1. 
826»j, 6 -0.7427<T 3 + 0.7427<r 4 ( 5 - 7 - 2 ) 

i) 3 = 0.24^ + 0.2949r; 2 - 0.9341 tj 3 + 0.2949r; 4 + 0.2949^ - 0.5901 »j 13 + 1.896»? 14 - 0.5901»j 15 

- 0.5901»j 16 - 0.24<t 3 + 0.24<t 4 - d, (5- 7 - 3 ) 

j) 4 = 0.001201tjj + 0.001475 tj 2 + 0.001475tj 3 — 1.228 tj 4 + 0.001475*7 5 — 0.002953»; 13 + 0.002953t/ 14 

+ 2.456tj 15 — 0.5901r; 1 g — 0.0012<7 3 + 0.0012c 4 + d 2 (5.7.4) 
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Table 1. Assignment of the Primitive (Perturbation) Variables to the Fundamental Variables 


6K 7 = 

u i» 

**21 = 

= u 8 , 

**22 = 

= u 9 , 

**27 = 

U 107 

**28 = 

: U Jl7 

* n l = 7?j, 

* n 2,2 = 

= di, 

* n 2,3 

= d 2 , 

* n 2,4 : 

- d 3 , 

*n 3 - ^37 

* n 9 ~ 

^12’ 

* n 10,2 = U 2’ 

6n 15 = 

: 

* n 17,l 

- u 4 , 

* n 17,2 

= u 5 , 

* n 17,3 


* n 17,4 

= U 7 , 

* n is = <r 4 , 

* n 27 = 

cr 17' 

* N 1,1 

- v 2 . 

™1,2 

= ^3’ 

* N 1,3 = 

--r] 4 . 

* N I,4 = 

- »?5> 

* N 5,i - 17> 

* N 5,2 : 

= 9s ’ 

* N 5,3 

= ^91 

* N 5,4 

= *?10’ 

6N 36,i 

= 9i3i 

^36,2 

= *?14> 

6N 36,3 - ^15 

™36,4 

= ^16’ 

* P 1 = 

<r 2’ 

* P 5 = 


6 p !2 = 


* P 13 = 

: a S' 

* P 14 - °9> 

* P 24 = 

: a 14> 

6P 25 : 

= ‘T'le* 

^ P 35 d 

= ff l. 

* p 36 - 

a lV 

* T 1 = 


*T 2 — d 4’ 

6T 4- 

^5’ 

* T 5 = 

^iit 

6T 7 = 

d 6 , 

* T 16 = 

^ 10 ’ 

«T 2 4 = 

- ^15' 

* T 35 = ff 6’ 

^ T 36 - 

: ^17’ 

iT 38 : 

= d 8 , 

* r 36 = 

°’l 2 > 

6t 37 = 

u 3’ 

V g = d 

I 7 . 

*V, = t? 19 , 


* w 37 - 


t?5 = 0.0161T7J + 0.01978»? 2 + 0.01978t? 3 + 0.1978r? 4 - 1.209 tj 5 - 0.03958t? 13 - 0.03958tj 14 

- 0.03958t? 15 + 2.509r? 16 - 0.0161<r 3 - 0.9839<r 4 + d 3 (5.7.5) 

1731.0^ 6 = 540.127?! - 7121. 0r? 6 - 583.4 <t 3 + 2981<t 6 - 648.8dj - 741d 2 - 742.4d 3 + 4080.0d 4 
+ 60.0d 5 — 133.3d 6 (5.7.6) 

2.0<r 2 = 2437. 0 t? 2 + 2437.0r? 3 + 2437.0 t? 4 + 2437.0t? 5 + 692.1.0tj 6 - 101,460.08r? 19 (5.7.7) 

«r 2 - .08688 <t 3 - <r 5 = 0 (5.7.8) 

i} 7 = 0.3163t? 2 - 0.9127t? 3 - 0.9127t? 4 - 0.9127t? 5 - 0.93»? 7 + 2.684i? 8 + 2.684t? 9 + 2.684r? 10 

- 0.7427t? 12 + 0.7427<r 3 - 0.7427u 2 (5.7.9) 


f ] 8 = - 0.2949r? 2 + 0.934 1?? 3 - 0.2949r? 4 - 0.2949 t? 5 + 0.8673t? 7 - 2.747ij 8 + 0.8673r? 9 + 0.8673r/, o 


- 0.24t? 12 + 0.24<r 3 - 0.24u 2 


(5.7.10) 


i} 9 = - 0.001475t? 2 - 0.001475t? 3 + 1.228r? 4 - 0.001475 t? 5 + 0.004324t? 7 + 0.004324t? 8 - 3.610t 7 9 
+ 0.004324r? lo - 0.0012 tj 12 + 0.0012<r 3 - 0.0012u 2 (5.7.1 1 ) 
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i) 10 = - 0.01978^2 - 0.01978t/ 3 - 0.01978»/ 4 + 3.556tj 5 + 0.05818»/ 7 + 0.05818r/ g + 0.05818 t/ 9 
- 3.556»/ 10 - 0.0161t; 12 + 0.0161<t 3 - 0.0161u 2 (5.7.12) 

Table 2. Assignment of the Fundamental Variables to the Primitive (Perturbation) Variables. 

dj = 6n 2 2 , d 2 = 6n 2 3 , d 3 = 6n 2 4 , d 4 = «5T 2 , d 5 - 6T 4' d « ~ 6T 7’ 

d 7 = \ g , d 8 = 6T 38 , r?,= 6 n v r> 2 = W lv ^“ 6N 1,3* 

*? 5 = ^Ni 4 , i/ 6 = 6Tj, »> 7 = 6N 5 i, '?8 = ^ n 5,2’ , ?9 = 5N 5,3’ ’ho - ^ N 5,4’ 

»>n- iT 5- »?i2- 6n 9- »7i3-^ N 36 t i’ ^4“ 5N 36.2’ ^5-^ N 36.3’ ^16- 6N 36, 4’ 

Vl7 = «T 36 , t , 18 = 6 w 37 , ^9 = ^, ^1=^35. <r 2 ±SP v <r 3 = 6n 3 , 

*16-* P 25’ ff 4~ in i8. <hs-* P 5’ ff 6- iT 35’ *7-^12’ ^8-^13’ 

<t 9 = 6P 14 , trio = ^ T 16’ a n “ ^ P 36’ ^12-^36’ a i3~ 6n ^ "u - SP 24' 

<7 15 = 6T 2 4, «r 17 = 6n 27 , u } = 6K 7 , u 2 = to 102 , u 3 = 6t 37 , u 4 = 6n 171 , 

u 5 = 5n 17 2 , u 6 = 6n 17 3 , = 6n 17 4 , u 8 = 6K 21 , u 9 = 6A 22 , u io — 6A 27 , 

Un = ^^28’ 

588.6^ jj = - 2981t/ 6 - 3041.0i/ n - 582.4 tj 12 + 583.4<r 3 - 582.4 u 2 + 60.0d 7 (5.7.13) 

i 12 = 1.279 <r 8 - 1.279 <r 9 (5.7.14) 

tr 5 - <t 7 = 35.9t ? 12 + 10,470 Ul + 35.9u 2 (5.7.15) 

. ... (5.7.16) 

<r 7 -<r 8 = 4.145r? 12 v 

0.68<t 5 = 2437.0i/ 7 + 2437.0»/ 8 + 2437.0t/ 9 + 2437 .0 »/ 10 + 235.4i/ n (5.7.17) 

0.02402<r 9 + 29.14<r 10 - 0.02402<r n = 29.14r/u + 92.48»/ 12 (5.7.18) 

2.457<r 9 - 2.457<r n + 54.19 <t 12 = 92.48»/ 12 - 62.82t/ 18 (5.7.19) 

rj ls = - 0.4051j/ 18 - 0.1725 a n + 0.1725u 3 (5.7.20) 

- 20.29<t 4 + 20.29<r 13 + 102.3<r ls = 0.42»?i + 100.2r/ 17 (5.7.21) 
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a l4- ff i6= !58.8 t7i (5.7.22) 

— <Tj + o’jg + 51.31<Tj 7 = 51.31 t/j + 658.3u n (5.7.23) 

“ a \ + ^16 “ 20.88<r 17 = 5874. Ou 10 (5.7.24) 

102.3<r 6 - 25.66<t 15 + 9.15cr 17 = 6.87r?j (5.7.25) 

- 16.26<7 4 + a u + 16.26<t 13 - a 14 = 16.26r7j + 839.4u g (5.7.26) 

cr n -763.9a 13 -<r 14 = 4.385ug (5.7.27) 

<r 4 = 0.01064(^3 75.7.281 


0 . 99870 -jj _ 2481 . 0 »? i 3 + 2481 . 0 jj 13 + 2481 - 0 r 7 14 + 2481 . 0;? 15 + 2481 . 0;? 16 + 345 . 9;? 17 ( 5 . 7 . 29 ) 


^13 - ~ 0 . 7472 r ; 1 + 0 . 93;? 7 - 2 . 684;? 8 - 2 . 684 t ? 9 - 2 . 684;? 10 + 0 . 7427;? 12 - 0 . 622 q 13 + 1 . 837;? 14 
+ 1 . 837 » 7 15 + 1 . 837 t 7 16 - 0 . 7472 a 4 + 0 . 7427 u 2 + u 4 ( 5 . 7 . 30 ) 

»? 14 = - 0 . 247 ?! - 0 . 8673 t ? 7 + 2 . 747 t ? 8 - 0 . 8673 t ? 9 - 0 . 8673 ;? lo + 0 . 24 t ? 12 + 0 . 5901 t ? 13 - 1 . 869 t ? 14 
+ 0 . 5901 t ? 15 + 0 . 5901 t ? 16 + 0.24 a 4 — 0 . 76 u 2 + u 4 ( 5 . 7 . 31 ) 

»? 15 = - 0 . 0012 ;?! - 0 . 004336 ;?- - 0.004336t? 8 + 3 . 61 r ? 9 - 0 . 004336 ;? lo + 0 . 0012;? 12 + 0 . 002951 f ?, 3 
+ 0 . 002951 r ? 14 — 2.455t?j 5 + 0.002951t? 16 — 0.0012<t 4 + 0 . 0012 u 2 + u 6 ( 5 . 7 . 32 ) 


*?16 - “ 0.01617?! - 0.05818;?- - 0.05818*? 8 - 0.05818;? 9 + 3.556;? 10 + 0.0161;? 12 + 0.03958t? 13 
+ 0.03958^i 4 + 0.03958;? 15 - 2.418 ^j 6 - 0.0161<r 4 + 0.0161u 2 + u 7 I 


1213.0;? 17 = - 8696.0;?i + 8696.0 ^j 2 - 2981.0t? 17 - 8696.0<t 4 + 2981.0<r 10 + 8708.0u 4 + 8800.0u 5 
+ ll,060.0u e + ll,070.0u 7 + 0.08597d g ( 5 7 


a u ~ a 9 = 4 °. 6 ;? i 8 + 1 . 792 ^ i 2 


(5.7.35) 
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^19 — ^7 


(5.7.36) 


5.7 - THE STATE VARIABLE FORM OF THE LINEARIZED 
MATHEMATICAL MODEL 

The equations are arranged compatibly with the symbolic form of the 
primitive equations from section 5.1. To obtain the state variable representation of 
the model (Equations 5.1.5 and 5.1.6), one can load the coefficients in the 36 
equations from Section 5.7, (which are in the form of Equations 5.1.1 and 5.1.2) 
into arrays Al, A2, A3, A4, A5, Bl, B2, B3, and B4 in a MATLAB compatible 
format for reduction to state space format. A compact MATLAB input form for 
each of the matrices follows: 

Al 

//coefficient of d(eta)/dt for study condition 1.; 

//25J189; 

al=eye(19); 

al(6,6)=1731; 

al(ll,ll)=588.6; 

al(17,17)=1213.; 

A2 

//A2 for study condition 1.; 

//24J189; 

a2(19,19)=0.; 

a2(2,l:5)=<.7427, -.3163, .9127,. 9127, .9127>; 
a2(2,13:16)=<.633, -1.826, -1.826, -1.826>; 
a2(3,l:5)=<.24, .2949, -.9341,. 2949, .2949>; 
a2(3,13:16)=<-.5901, 1.869, -.5901, -.5901>; 

a2(4,l:5)=<.001201, .001475,-001475, -1.228, .001475>; 
a2(4,13:16)=<-.002953, -.002953, 2.456, -.002953>; 
a2(5,l:5)=<.0161, .01978,-01978,. 01978, -1.209>; 
a2(5,13:16)=<-.03958, -.03958, -.03958, 2.509>; 
a2(6,l)=540.1;a2(6,6)=-7121.; 

a2(7,2:12)=<.3163, -.9127, -.9127, -.9127,0., -.93, 2.684,2.684, 2.684,... 

0,-.7427>; 

a2(8,2:12)=<-.2949, .9341, -.2949, -.2949,0.,. 8673, -2.747,... 

.8673, . 8693, 0.,-.24>; 

a2(9,2:12)=<-.001475, -.001475, 1.228, -.001475,0., .004324, .004324,... 

-3.61,-004324, 0.,-.0012>; 

a2(10,2:12)=<-.01978, -.01978, -.01978, 1.209,0., .05818, .05818,. 05818,... 

-3.556, 0.,-.0161>; 
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a2(ll,6)=2981.;a2(ll,ll)=-3041.;a2(ll,12)=-582.4; 

a2(13,l)=-.7472;a2(13,7:16)=<.93, -2.68, -2.684, -2.684,0.,... 

.7427, -.622, 1.837, 1.837, 1.837>; 
a2(14,l)=-.24;a2( 14,7: 16)= <-.8673, 2. 747,-. 8673,-. 8673,0.,... 

.24, .5901, -1.869, .5901,. 5901>; 

a2( 15,1)=-. 00 12;a2( 15,7: 16)= <-.004336, -.004336,3. 61, -.004336,0.,... 

.0012, .002951, .002951, -2.455,. 002951>; 
a2(16,l)=-.0161;a2(16,7:16)=<-.05818, -.05818, -.05818,3.556,... 

0.,.0161,. 03958, .03958,. 03958, -2.418>; 
a2(17,l)=-8696.;a2(17,12)=8676.;a2(17,17)=-2981.; 
a2(18, 18)=-. 04051; 

A3 

//This is A3 for control study model 1. 

a3(19,17)=0., 

a3(l,l:2)=<1.279,-1.279>; 

a3(2,3:4)=<-.7427,.7427>; 

a3(3,3:4)=<-.24,.24>; 

a3(4,3:4)=<-.0012,.0012>; 

a3(5,3:4)=<-.0161,-.9839>; 

a3(6,3)=-583.4;a3(6,6)=2981.; 

a3(7,3)=.7427; 

a3(8,3)=.24; 

a3(9,3)=.0012; 

a3(10,3)=.0161; 

a3(ll,3)=-583.4; 

a3(12,8:9)=<1.279,-1.279>; 

a3(13,4)=-.7424; 

a3(14,4)=.24; 

a3(15,4)=-.0012; 

a3(16,4)=-.0161; 

a3(17,4)=-8696.;a3(17,10)=2981.; 

a3(18,12)=-.1725; 

A4 

//This is A4 for control study case 1; 

//24J189; 

a4(19,ll)=0.; 

a4(7,2)=-.7427; 

a4(8,2)=-.24; 

a4(9,2)=-.0012; 

a4(10,2)=-.0161; 

a4(ll,2)=-582.4; 

a4(13,2:4)=<.7427,0.,l.>; 

a4(14,2)=-.76;a4(14,5)=l.; 

a4( 15,2)=. 001 2;a4(15,6)=l.; 

a4(16,2)=.0161;a4(16,7)=l.; 

a4( 17, 4:7)=<8708., 8800., 11 060., 10070. >; 

a4(18,3)=.1725; 

A5 

//This is A5 for control study system 1.; 
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//24J189; 

a5(19,8)=0.; 

a5(3,l)=-l; 

a5(4,2)=-l.; 

a5(5,3)=l.; 

a5(6!l:6)=’<-648.8, -741, -742.1, 4080., 60., -133.3>; 

a5(ll,7)=60.; 

a5(17,8)=. 08597; 

a5(19,6)=l.; 

B1 

//this is the B1 matrix for the first control study model; 

//24J189; 

bl(17,17)=0.; 

bl(l,2)=2.; 

bl(2,2)=l.;bl(2,3)=-8.688e-2;bl(2,5)=-l; 

bl(3,5)=l.;bl(3,7)=-l; 

bl(4,7)=l;bl(4,8)=-l; 

bl(5,5)=.68; 

bl(6,9)=.02402;bl(6,10)=29.14;bl(6,ll)=-.02402; 

bl(7,9)=2.457;bl(7,ll)=-2.457;bl(7,12)=54.19; 

bl(8,ll)=.9987; 

bl(9,4)=-16.26;bl(9,ll)=l.;bl(9,13)=16.26;bl(9,14)=-l.; 

bl(10,ll)=l.;bl(10,13)=-763.9;bl(10,14)=-l.; 

bl(ll,4)=-20.29;bl(ll,13)=20.29;bl(ll,15)=102.3; 

bl(12,14)=l.;bl(12,16)=-l.; 

bl(13,l)=-l.;bl(13,16)=l.;bl(13,17)=51.31; 

bl(14,6)=102.3;bl(14,15)=-76.64;bl(14,17)=9.; 

bl(15,4)=l.;bl(15,13)=-.01064; 

bl(16,9)=-l.;bl(16,ll)=l.; 

bl(17,l)=-l.;bl(17,16)=l.;bl(17,17)=-20.88; 


B2 

//this is b2 for the first control studies model.; 

//24JI89; 

b2(17,19)=0.; 

b2(l, 2:6)=<2437., 2437. ,2437. ,2437., 692. I>;b2(l,19)=-101460.; 

b2(3,12)=35.9; 

b2(4,12)=4.145; 

b2(5,7:ll)=<2437, 2437, 2437, 2437, 235.4>; 
b2(6,ll)=29.14, 

b2(7,12)=92.48;b2(7,18)=-62.82; 
b2(8,13:17)=<2481, 2481, 2481, 2481, 345.9>; 
b2(9,l)=16.26; 

b2(ll,l)=.42;b2(ll,17)=100.2; 

b2(12,l)=158.8; 

b2(13,l)=51.31; 

b2(14,l)=-6.87; 

b2(16,12)=1.792;b2(16,18)=40.60; 

B3 

//this is b3 for control studies model 1; 
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//24J189; 

b3(17,ll)=0.; 

b3(3,l:2)=< 10470, 35.9>; 

b3(9,8)=839.4; 

b3(10,9)=4.385; 

b3(13,ll)=658.3; 

b3(17,10)=5874.; 

B4 

//B4 is empty 

6 - SOME PROPERTIES OF THE LINEARIZED STATE VARIABLE MODEL 

With the primitive form matrices available for manipulation via MATLAB, 
the determination of the linear state variable representation coefficient matrices is 
quite straightforward. SVD analysis shows both A, and Bj to be maximum rank, 
and the resulting A, % C, % 8, and <J matrices in Equations 5.1.5 and 5.1.6 are 
quickly generated. For a system of this size, the specific numerical values of any 
element is of tertiary significance, and thus the presentation of these matrices is 
relegated to Appendix 3. However, an item of great interest is the eigenvalues of 
the state coefficient matrix, A: 


EVV = 

1.0D+04 * 
-5.266821895634872 
-0. 021360198075745 
-0. 003977263492918 
-0. 003977263492918 
-0. 000459990795721 
-0. 000459990795721 
-0. 000093892720314 
-0. 000021248217594 
-0. 000366233205798 
-0. 000365106096282 
-0. 000365101237584 
-0. 000366233205798 
-0. 000365106096282 
-0. 000365101237584 
0.000002445599342 
0.000000294968294 
-0. 000000000002605 
0.000000027480837 


O.OOOOOOOOOOOOOOOi 
O.OOOOOOOOOOOOOOOi 
0.009552095834877i 
-0.009552095834877i 
0.000120566056909i 
-0.000120566056909i 
-O.OOOOOOOOOOOOOOOi 
-O.OOOOOOOOOOOOOOOi 
0.000176262166536i 
0.000173786717458i 
0.000173829082129i 
-0.0001 762621665361 
-0.000173786717458i 
-0.000 173829082 129i 
O.OOOOOOOOOOOOOOOi 
O.OOOOOOOOOOOOOOOi 
O.OOOOOOOOOOOOOOOi 
O.OOOOOOOOOOOOOOOi 
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EW(continued) 

0.000000000000000 O.OOOOOOOOOOOOOOOi 

It is seen that the four left-most eigenvalues are greater than 39.0 in 
magnitude, suggesting that a representation of 15th order (a reduced order 
representation) should be acceptable for control studies. The remainder reflect 
dynamics which will require at least 2 seconds to play out, and which may be very 
important to the determination of the feedback control structure. One should be 
mindful that the very small (but non-zero) eigenvalues might be due to numerical 
error in the computational finite arithmetic, but a singular value analysis indicates 
that in this case the small but non-zero eigenvalues are in fact just that, and the 
plant is very mildly (slow dynamics) unstable about the prescribed experimental 
operating condition. This conclusion is consistent with the physical nature of the 
system. To further understand the dynamic behavior of the model, it is important 
to reformulate the model into smaller individual subsystems, the canonical basis. 
We employ the Real Jordan canonical decomposition. 


6.1 - A REAL JORDAN CANONICAL FORM OF THE LINEAR STATE 
VARIABLE REPRESENTATION 

For systems of this dimension, the Jordan canonical decomposition is of 
great value to gain insights into the system behavior (see Brogan, 1985). In the 
classical Complex Jordan canonical transformation (JCT), one can investigate the 
system in the canonical state , which is characterized by 

*(t) = MJq(t) 

6j(t) = Aq C j(t) + BJ*(t) + CJHt) 

where 


54 



MJ is the Complex Jordan modal matrix (CJMM) (generally complex) 

9 cj{t) Complex Jordan state vector (generally complex) 

A is the Complex Jordan eigenvalue matrix (CJEM) (generally complex and 
sometimes 

diagonal but often having some ones in the first super diagonal) 

BJ = MJ ~ 1< 3B (generally complex), and 
CJ = MJ ~ l C (generally complex). 

Often the JCT is satisfactory, but in the present case the JCT produces a 
CJMM requiring thirty-eight (38) columns to present, half of them imaginary. By 
employing a real Jordan Canonical transformation (RJCT) (see Takahashi, 
Rabins, and Auslander, 1970), one can obtain a nineteen column Real Jordan 
Modal Matrix (RJMM) MR and a Real Jordan Eigenvalue Matrix (RJEM) ER. If 
any entry in A is complex, then the Real Jordan canonical state vector differs from 
the Complex Jordan state vector, although they axe related. In the sequel we 
employ the asymmetric Real Jordan Eigenvalue Matrix (RJEM) form (Takahashi, 
et al, 1970). To denote this case, 


*(<) = MRq(t) 

«(<) = ERat) + BRu(t) + CRJ(t) 

where 

MR is the Real Jordan modal matrix (real) 

ER is the eigenvalue matrix (real and tri-diagonal possibly with block 
asymmetry; in some cases having some ones on the superdiagonal just above 
the diagonal) 

BR = MR~ 1( $ (real), and 
CR = MR-'C (real). 

Four forms of the Real Jordan Transformation are available (Takahashi, et 
al i 1970, and Siljak, 1986. In the version we use here Takahashi, et al, 1970), the 
resulting Real Jordan Eigenvalue Matrix is tri-diagonal and block asymmetric. The 
real part of each of the complex Jordan eigenvalues appears on the diagonal, and 


55 



the negative imaginary part appears above the diagonal while its negative appears 
below the diagonal. Here, to save space, we present the blocks concatenated in a 
two-row, nineteen column matrix; 


DBER = 

1.0D+04 * 

Columns 1 thru 4 

-5.266821895634871 0.000000000000000 

0.000000000000000 -0.021360198075745 
Columns 5 thru 8 

-0.000459990795721 -0.000120566056908 
0.000120566056908 -0.000459990795721 
Columns 9 thru 12 

-0.000366233205798 0.000176262166536 

-0.000176262166536 -0.000366233205798 
Columns 13 thru 16 

-0.000365 1 01237584 0 .000 173829082 129 
-0.000173829082129 -0.000365101237584 
Columns 17 thru 19 

-0.000000000002605 0.000000000000000 

0.000000000000000 0.000000027480838 


-0.003977263492919 0.009552095834876 
-0.009552095834876 -0.0039772634929 1 8 

-0.000093892720314 0.000000000000000 

0.000000000000000 -0.000021248217594 

-0.000365106096282 0.000173786717458 

-0.000173786717458 -0.000365106096282 

0.000002445599342 0.000000000000000 
0.000000000000000 0.000000294968294 

0.000000000000000 

0.000000000000000 


The associated Read Jordan modal matrix is 


MR = 

Columns 1 thru 4 

0.021374485835695 0.98457631 1054076 
0.487360549065908 -0.00 1955648 123576 
0.157488261878300 -0.000631875128675 
0.000787441309857 -0.000003163992840 
0.010564837649393 -0.00004 1386247457 
0.221 177235667850 0.000096446984620 
-0.487360123124829 -0.000868067795239 
-0.157488125006523 -0.000280512407754 
-0.000787440784020 -0.00000 1402604758 
-0.010564828397520 -0.000018823588278 
0.650447775541264 0.002825999155352 

0.051688720402714 0.172709078715878 

-0.0000004 12 15397 1 0.00284426633072 1 

-0.000000 136513800 0.000911896962854 
-0.000000000675574 0.000004561937004 

-0.000000009 161841 0.00006 1206 155787 

-0.000034465044499 0.027550330062014 

0.000000302667457 0.000249544898 103 

0.000000000000000 0.000000000000000 
Columns 5 thru 8 

0.023978462944404 0.029493847479883 

-0.007084457634296 0.02531 1329129724 

-0.002320135705025 0.008169161137254 


-0.222861859233880 0.322197951244839 
0.005816754894363 -0.003751480241953 
0.00 1 879634406395 -0.0012 12404055 1 58 
0.000009401663918 -0.000006061316275 
0.000125394594452 -0.000081420370778 
0.002584220507461 -0.001201962066472 
0.001402104579983 -0.000485856568940 
0.000453080327699 -0.000157002400 1 03 
0.000002265548592 -0.000000784867518 
0.000030396185574 -0.000010523959557 
0.010702876866564 -0.012052105651288 
-0.381326891771549 -0.833086493557783 
-0.007234575766325 0.004234483907949 

-0.002332320942 123 0.00 1 369490254 1 80 
-0.00001 1663425357 0.000006847095427 
-0.000156484554555 0.000091864777795 
-0.071485556002443 0.039486759722147 
0.00 1 859344897543 -0.002002273880364 
0.000000000000000 0.000000000000000 

-0.002568503774483 - 0 . 1703336237 13372 
0.016182368235621 0.009271856159084 

0.005169242141229 0.003569791470073 
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MR(continued) 
-0.000012098438914 
-0.000164258756038 
0.202604293203552 
-0.025478334763020 
-0.008215011499415 
-0.000041406652986 
-0.000538211083104 
0.523635785674718 
-0.138683675631643 
0.032485032408920 
0.010534793286127 
0.000053380935448 
0.000688589834673 
-0.132208556882285 
-0.011987363542998 
0.000000000000000 
Columns 9 thru 
-0.000434549338021 
-0.333604970041021 
-0.110676221670387 
-0.000539486289905 
0.450460633757361 
0.004695728699299 
0.069447849286962 
0.011517247536669 
0.000113388994367 
-0.076749231164386 
-0.020268951933378 
0.002413053767813 
0.264173290919339 
0.099126440708361 
0.000425958175786 
-0.362464456254652 
0.015262522898163 
-0.000309454575899 
0.000000000000000 
Columns 13 thru 
-0.000078989144497 
0.009929433211928 
0.014324858119588 
-0.024252456132820 
0.000038235837291 
-0.000397613710317 
-0.264356820743857 
-0.045696087138582 
0.310129833762039 
0.000019931422614 
-0.001259213077060 
0.000452642963183 
0.254464536530617 
0.031371215925709 


0.000040819503345 

0.000533867262298 

-0.038048117289116 

-0.049290720473403 

-0.015943111731417 

-0.000080134199729 

-0.001067719061769 

0.769014585609127 

-0.172002174168538 

0.023965232663896 

0.007773824942749 

0.000039173375510 

0.000528457533651 

-0.134156712888075 

-0.008688571309677 

0.000000000000000 

12 

-0.002365889272441 

0.083748057550945 

0.011827578249793 

0.000137051277497 

-0.098016095540404 

-0.012815136031644 

0.234874927207185 

0.078290731338212 

0.000379988557379 

-0.311843457165495 

-0.038544150768501 

0.013556551929108 

-0.318691699740971 

-0.089963600874162 

-0.000516962540559 

0.405204053935332 

0.006277594254333 

0.001037430467572 

0.000000000000000 

16 

-0.000001534446000 

-0.373765232568121 

-0.064326976103987 

0.438160068345558 

0.000031136207611 

-0.000268080980026 

-0.001438689585612 

-0.009178093458171 

0.010607832143095 

-0.000028474169719 

0.000467889438399 

0.000012798653809 

0.375231370860600 

0.073505078970693 


0.000025905961077 

0.000362212241106 

0.171601361074757 

0.003741219057719 

0.001181558255228 

0.000005911968398 

0.000086347194217 

0.196136838947666 

0.022504076070899 

-0.019631368541019 

-0.006350693074713 

-0.000031773551538 

-0.000409479782940 

0.436058794716317 

0.008888177603467 

0.000000000000000 

-0.000168637166761 

-0.324410007385080 

0.323082733665237 

0.001166063755237 

0.000216340950950 

-0.000845902687344 

0.183247014243786 

-0.182999358919719 

0.000033286646231 

-0.000080476569712 

-0.002726952896061 

0.000966820438091 

0.141231439230398 

-0.140083367030527 

-0.001199002422559 

-0.000130725847818 

0.000560853205478 

0.000066645761575 

0.000000000000000 

0.002127499073159 

-0.118764556142297 

-0.013039260699884 

0.000108773981237 

-0.850335665841901 

-0.004624688434761 

-0.041878623768547 

-0.005138670800554 

0.000035151784017 

-0.287200647192445 

-0.001593536101620 

-0.003764851461600 

-0.056307401242189 

-0.005299944868085 


0.000017923850718 

0.000196493162171 

0.040494054146723 

-0.001732833228306 

-0.000352542308634 

-0.000001727130129 

-0.000039178193928 

0.108489631181109 

-0.168037356184343 

-0.010424941737283 

-0.003131532280652 

-0.000015571595920 

-0.000234635785042 

0.092275627607291 

-0.952601895329476 

0.000000000000000 

-0.000007487154907 

0.252277615693271 

-0.252065311662433 

0.000071652116033 

-0.000110198673781' 

-0.000564650385160 

0.225511990573353 

-0.224609147353368 

-0.000825646783815 

-0.000151242264971 

0.000813576716260 

0.000050523205025 

-0.477511486074352 

0.476674545754798 

0.000753955953713 

0.000257469144136 

-0.001209068228379 

0.000037619146834 

0.000000000000000 

-0.001249937877648 

-0.857448582488209 

-0.001038607482383 

0.000877705234494 

-0.000277775049492 

0.003685286338719 

-0.291430937702131 

-0.000385007648221 

0.000299289112009 

-0.000097086177615 

0.002980417496428 

0.002041513814356 

-0.423848960773323 

0.001312670859152 
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MR( continued) 
-0.285862028060663 
-0.000057032339740 
0.000298663976310 
0.000031581325941 
0.000000000000000 


-0.448637501247626 

-0.000002070183515 

-0.000615576198059 

0.000016846781909 

0.000000000000000 


0.000061912601134 

-0.418055811773158 

-0.018368371512324 

0.006360851231764 

0.000000000000000 


0.000448816238949 

-0.000015706546943 

0.012426412413986 

-0.003909868113154 

0.000000000000000 


Columns 17 thru 
0.000016412362758 
-0.668318523019703 
0.534857523446167 
-0.001354064207911 
-0.000041012067976 
-0.000059372032713 
-0.227253876620181 
0.181877033253038 
-0.000460405901928 
-0.000013694027642 
-0.000060628271828 
-0.000026516334553 
-0.333067189232513 
0.267615210333714 
-0.000675540588126 
0.000000040019953 
-0.000183998003297 
0.000051731173235 
0.000000000000000 


19 

0.000023656694458 

-0.634070002293218 

-0.000070529334320 

0.574758903861677 

-0.000014753075049 

-0.000023887458395 

-0.215654455358315 

-0.000038953357662 

0.195499349292086 

-0.000006059271029 

0.000031561837089 

-0.000038039401422 

-0.316633187948755 

0.000074545807038 

0.287443866902551 

-0.000000005204576 

-0.000145735125755 

0.000074082924124 

0.000000000000000 


0.000000051400872 

0.950909182004328 

0.307243950675371 

0.001536156292617 

0.020611724800414 

0.000036634532086 

0.000000533624518 

-0.000000411699438 

0.000000001131302 

0.000000000197389 

0.000070776809305 

0.000000051505625 

0.000000776667388 

-0.000000607552454 

0.000000001651331 

0.000000000125281 

0.000070773448986 

-0.000000100483148 

0.030751715660749 


The rows of MR show how the Jordan canonical state variables contribute 
to the descriptive state variables. As an example, it is clear that Jordan state 19 
solely determines descriptive state 19 (the volume of the upper chamber). All the 
others are not so clear-cut, and this is generally the case. Another item of interest 
is the influence the individual processes of the primitive form have on the Jordan 


elements, and how the descriptive state variables contribute to the Jordan 


canonical states. 


The first question is not so easily answered, and strictly involves (1) 
perturbing the coefficients of the individual process equations one at a time, (2) 
constructing the A it A 2 , ■ • • and B v B v • ■ • , (3) constructing .X, and (4) obtaining 
the Jordan elements, for each perturbation ; a very large task which in some cases may 
be well worth the effort. 
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The second question is much more easily addressed; one needs the effective 

inverse of MR, MRNV. 

MRNV = 

1.0D+02 * 

Columns 1 thru 4 

-0.000000079360991 0.004058094461690 0.004058094461691 0.004058094605937 
0.011027036117855 0.046300274368076 0.046300258010397 0.046300248492862 
0.004303637906307 0.133032132081097 0.133032326987938 0.133032162106799 
0.000330492422439 -0.047947626278134 -0.047947545627571 -0.047947473501417 
-0.000028472726969 -0.021029860724648 -0.020986976831392 -0.021038110671024 
-0.000064132701893 0.003434273108404 0.003446348835286 0.003424352323812 
-0.000051803236422 0.017930573774850 0.017974005745429 0.017960546499401 
0.000011053728587 0.000799112951901 0.000806778357354 0.000803074701879 
-0.000000014223569 -0.000174749685292 -0.000181435760722 -0.000176206203828 
-0.000000001495674 -0.000051705146068 -0.000052420806710 -0.000052802519075 
-0.000000029038514 -0.003065888425097 0.009491348409376 -0.001254296004910 
0.000000062645297 0.000496610685218 -0.001505837137098 -0.000107602331880 
-0.000000000083374 0.000009059687593 0.000003534084696 -0.006296304394037 
-0.000000000178431 -0.000004170475201 -0.000037576131655 0.010064370948144 
0.000000003635636 0.000001036616354 0.000001080948813 0.000001074283795 
0.000000069891068 -0.006362957631631 -0.007960180072798 -0.007016244971145 
0.000000002864779 -0.000000794483606 0.010158017373222 -0.000000532029538 
-0.000000000084737 0.000009776455899 0.000036147127610 0.009464088565077 
0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 
Columns 5 thru 8 

0.004058094394871 0.001152686400887 -0.011935557140388 -0.011935557208376 
0.046300285905894 0.013265246187142 0.057116207072712 0.057119033703496 
0.133032135860233 0.037911355587235 0.160928666851032 0.160928822588576 
-0.047947738407967 -0.014374467244194 -0.061640535666758 -0.061658213660490 
-0.020691099176059 0.034921422890468 -0.011098488099617 -0.011227956679605 
0.003803974812573 -0.027779428033610 0.007250430113147 0.007259378246052 
0.017254737184331 0.008719287158453 0.027215376969202 0.027225627932003 
0.000712572446781 0.000002971673030 0.000813120253453 0.000819460384678 
0.010710862653134 0.000000005468609 0.000270145171574 0.000282356727909 
0.003122887569759 0.000001590921333 0.000349617405094 0.000360623476204 
-0.000063017657065 0.000003234124340 0.006737159169779 -0.020841670347579 
-0.000314719087115 0.000000103145409 0.003360182233558 -0.010450036641585 
-0.000001308118171 -0.000000011508949 -0.000023667284594 -0.000060026655955 
-0.000000122913299 0.000000038795494 -0.000004712455155 0.000068936054341 
-0.006360029205293 -0.000000008741358 0.000003202670888 0.000003247885522 
0.000998588949181 -0.000000040047863 -0.006361272919333 -0.007962328802177 
-0.000151815240640 0.000000000005285 -0.000000794104094 0.010158017537432 
-0.000000656698796 0.000000000012694 0.000009772131394 0.000036148699686 
0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 
Columns 9 thru 12 

-0.011935557213817 -0.011935556727441 -0.001153018373085 -0.000000248506030 
0.057116497941355 0.057116149439601 0.005877579478471 0.002019520710580 
0.160928817515437 0.160924826314253 0.015481543135955 -0.004711091034897 
-0.061642367708082 -0.061639763771372 -0.008191542717668 -0.009418427332986 
-0.011069980331601 -0.012388839824255 -0.001926214702688 0.000040158286908 
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MRNV(continued) 

0.007268062053405 0.006569233515909 

0.027218872955276 0.027057618485695 

0.000816391908286 0.000741546948676 

0.000271282179571 -0.016734526561433 
0.000353721933110 -0.021552321505284 
0.002311184512667 -0.000328079037130 
0.00 1 98844760 1 7 1 6 0.0007 1 309 1 7 1 3345 

0.026630590408530 0.0000024902656 15 
-0.010920585174335 0.000001975484830 

0.000003240250745 -0.006484464908057 
-0.007016128237641 0.001017945335061 

-0.000000531650199 -0.000151814857910 
0.009466197440280 -0.000000663357369 
0.000000000000000 0.000000000000000 
Columns 13 thru 16 

-0.0000000 10232328 -0.0000000 1 0240033 
-0.133980202090320 -0.133980287282974 
-0.38277 1 1 26779035 -0.38277 1142181409 
0.140439624292678 0.140440163283337 

0.05045 143 1 376845 0.050422807895227 

-0.012030447082541 -0.012061 160764686 
-0.054798544384965 -0.054723502606230 
-0.002087890892367 -0.002079371998441 
0.000165692758078 0.000170723208760 

-0.000 133847727 1 89 -0.000 1 39761062429 
0.001550864973431 -0.004810297437559 
-0.003270455776704 0.010134931293734 

-0.000001985481366 0.000033805350703 

0.000011478775925 0.000028203612623 

0.00000 1608382963 0.00000 165225585 1 
-0.006372995381737 -0.007968303847082 
-0.000000373998128 0.010158437949212 

0.00000980 1360430 0.000036 16908449 1 

0.000000000000000 0.000000000000000 


0.0 13288497480227 -0.000100749 15 1063 
0.006410068997662 -0.000068165259628 
0.00001 1816052571 0.000010907601423 

0.000000402932000 -0.000000012114450 
-0.000000127640973 0.0000000021 18973 
-0.000000279987031 -0.00000002 1 340049 
-0.0000008 1 8879 1 14 0.0000000589 11711 
0.000000016567716 -0.000000000077779 
-0.000000000840476 -0.000000000095633 
0.000000131861654 0.000000001504711 

0.000001063842386 0.000000054453148 

0.000000000007851 0.000000002793429 

-0.000000000410084 -0.000000000055846 
0.000000000000000 0.000000000000000 

-0.000000010228713 -0.000000017134799 
-0.133979525845905 -0.133998944505850 
-0.382771042151076 -0.382774278702508 
0. 140435389 120294 0. 14055746103783 1 
0.050432091026318 0.051033750070013 
-0.012023977207041 -0.012278165752610 
-0.054788422679860 -0.0549263002 17334 
-0.002084544597349 -0.002155565300914 
0.000167688144985 -0.010289394705628 
-0.000 1 344 1 295313 1 0.008450891278060 

0.000936745540861 0.000361908678533 

-0.001 127783648698 0.0001 599672779 1 8 
-0.005522465530611 0.000000996741591 

-0.012696905624139 -0.000001209177829 
0.000001646354553 -0.006529960873376 
-0.007028360747631 0.001016652104691 

-0.000000111581489 -0.000157015973551 
0.009466914237013 -0.000000654927085 
0.000000000000000 0.000000000000000 


Columns 17 thru 19 

-0.000000029080641 0.000000000245057 -0.168951276001464 

-0.018977198942561 -0.000484003880011 -1.927626665785087 
-0.053399192557579 0.005275456198033 -5.538549285503272 

0.021447381911251 -0.000419192722920 1.996212015962940 

-0.008633654390757 -0.000567608615953 0.874873693278677 

0.003955724022687 0.000950500410444 -0.143355254541355 

0 012028361877980 0.003266159255409 -0.746547286512787 

0.000030505602496 -0.010457639483299 -0.033288760591338 
0.000000086887107 0.000000133434329 0.000046090546549 
-0.000000466753684 -0.000000087561853 0.000032058549404 

-0 000000879334351 -0.000000076367255 0.000079066352446 

-0 000000214565499 -0.000000832746981 -0.000094870844419 
-0.000000000258151 0.000000001592361 -0.000000055535734 

-0.000000009969529 0.000000000864010 0.000001721555178 

0.000000375872476 0.000000597411480 0.004219980627829 

0.000002758842500 0.000018968559298 0.275968582994403 
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MRNV(continued) 

0.000000012863846 0.000000924297613 -0.101363240880237 
-0.000000005200491 -0.000000021013929 -0.001135783570108 
0.000000000000000 0.000000000000000 0.325185108704812 

The rows of MRNV show how the descriptive states map into the canonical 
states, e.g., a zero element shows that the corresponding descriptive state does not 
influence the canonical state. Again the graphic example of the present MRNV is 
that only descriptive state 19 influences canonical state 19. 


To a large degree, one can anticipate control problems by surveying the 


rows of BR. 


BR = 

1.0D+04 * 

Columns 1 thru 4 

0.000033277765579 0.000130878405621 
-0.270436643730725 -0.001556613041551 
0.630868324901464 0.000400677767752 
1.261233848115660 0.005022089592311 
-0.005377648405604 0.000 1 12434200427 
0.013491449792794 -0.000157447850720 
0.009128098531653 -0.000304870438139 
-0.001460651086450 -0.000013335794888 
0.000001622261810 -0.000001712921902 
-0.000000283754463 0.00000 143094 1 848 
0.000002857673731 0.000048166826014 
-0.00000788894386 1 -0.000 10 1336533629 
0.000000010415520 -0.000000338031286 
0.000000012806328 -0.000000281318062 
-0.000000201497800 -0.00000004 1532455 
-0.000007291891854 0.000079545751508 
-0.000000374071657 -0.000101582365427 
0.00000000747845 1 -0.00000036 1 385 1 78 
0.000000000000000 0.000000000000000 
Columns 5 thru 8 

-0.000000002212125 -0.000000002753828 
-0.002716549374846 -0.003070115293947 
-0.007701684171251 -0.008696589808871 
0.002960353494491 0.003359902480967 
-0.000122121126643 -0.000282886159496 
0.000166365897709 0.000240438774433 
0.000325391392126 0.000548848521515 
-0.0000 185806 1 7742 -0.0000 1 8063978837 
0.000001713535522 0.000001684803720 
-0.000001431472392 -0.000001386687617 
-0.000048166767799 0.000009287278674 
0.000101333746768 -0.000011297400332 


0.000000000000423 -0.000000002 1 89992 
-0.000000834906693 -0.002702155264034 
0.000009100161942 -0.007661 183393029 
-0.000000723107447 0.002944081335121 
-0.000000979124863 -0.000115286695586 
0.000001639613208 0.000163672815156 
0.000005634 1247 1 6 0.00031 55 1 806 1 785 
-0.000018039428109 -0.000018688943660 
0.000000000230 1 74 0.000001 663 165115 
-0. 000000000 15 1044 -0.00000 1 37 1985030 
-0.000000000131734 0.000015445523242 
-0.00000000 1 436489 -0.00003271 996 1 200 
0.000000000002747 -0.0000000 1 9873346 
0.000000000001490 0.000000114072057 
0.00000000 1 030535 0.00000004306732 1 
0.000000032720765 -0.00006353 1899403 
0.000000001594413 -0.000000002816499 
-0.000000000036249 0.000000097640266 
0.000000000000000 0.000000000000000 

-0.000000002585545 -0.000002306030737 
-0.002915425499070 -0.115956370297978 
-0.008260800240074 -0.045287281245884 
0.003 1 86078615706 -0.00346 1 762875487 
-0.000206405283429 0.000208279141379 
0.0002056 1191 9625 0.00074743546 1 862 
0.000449299274094 0.000514396152568 
-0.0000 1 9023159875 -0.000 116123010280 
-0.000102886733922 -0.000001067549952 
0.000084470 164144 -0.00000035392 1 73 1 
0.0000035460868 1 1 0.000000039440308 
0.00000 1 58 1860 1 28 -0.000000066 100955 
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BR( continued) 
0.000000338034779 
0.000000281312863 
0.000000043791131 
-0.000079482891612 
0.000101585312731 
0.000000361313563 
0.000000000000000 
Columns 9 thru 
0.000000016497528 
-0.000012689405687 
-0.000004786497825 
-0.000000451242909 
0.000000508824808 
-0.000000307751985 
0.000000218150669 
-0.000000013315908 
0.000000006374565 
0.000000001932760 
0.000000001422565 
-0.000000003168099 
-0.000000000010540 
-0.000000000007410 
-0.000000003808865 
0.000000002811293 
-0.000000003010574 
-0.000000000010687 
0.000000000000000 


-0.000055224678844 

-0.000126969965252 

0.000000050735183 

-0.000070032059265 

0.000000000057096 

0.000094668668195 

0.000000000000000 

11 

0.000146339651628 

-0.587192408382382 

-0.225134351923371 

-0.019419884191736 

0.005825476261093 

-0.000000023556181 

0.003841121645717 

-0.000589886727719 

0.000000760192096 

0.000000275993849 

0.000001949313747 

-0.000003332452360 

0.000000003033222 

0.000000014310614 

-0.000000195215325 

-0.000003737018697 

-0.000000152974373 

0.000000004526379 

0.000000000000000 


0.000000009945985 

-0.000000012919422 

-0.000065268404811 

0.000010395552737 

-0.000001569091814 

-0.000000006981002 

0.000000000000000 

-0.000015732112866 

-0.027037126660512 

-0.011004316774894 

-0.000606243705685 

-0.000413131342874 

0.000539977992367 

0.000005690443116 

-0.000026959915354 

0.000000034562735 

-0.000000018337603 

0.000000026034463 

-0.000000153983516 

0.000000000362044 

-0.000000000101466 

-0.000000008733001 

-0.000000169650612 

-0.000000006976594 

0.000000000206182 

0.000000000000000 


0.000000002912779 

0.000000003334174 

0.000000690081205 

-0.000001288497919 

0.000000545544026 

0.000000002955551 

0.000000000000000 


If the entries corresponding 


to a canonical block tire all 


zero, then that 


canonical state cannot be influenced by any of the control variables and it is said to 


be uncontrollable. If, in addition, the uncontrollable block is not stable, then all 
descriptive states which are influenced by the unstable canonical block cannot be 
controlled. If all entries in the set are not zero, but all have very small values, then 
one is faced with a condition of poor control authority, and the hardware 
realization of an effective controller system is likely to be difficult. 


CR = 

1.0D+02 * 

Columns 1 thru 4 

-0.004490135441989 -0.004551532285346 0.003563924216883 0.002716903821848 
-0.051272234744203 -0.051978785422193 0.040613319299609 0.031266438153402 

-0.147241967372107 -0.149261113285390 0.116779093121188 0.089357787865926 
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CR (continued) 
0.053335271998474 
0.007897999840441 
0.006965721013476 
-0.021242101359759 
-0.000807892176801 
0.000181433711021 
0.000051824509910 
-0.009492560598788 
0.001505798476935 
-0.000003529771001 
0.000037561590629 
-0.000001077672445 
0.007960195083228 
-0.010158017375203 
-0.000036147132368 
0.000000000000000 
0.000039954467968 
0.000459800561079 
0.001314085115675 
-0.000498248431341 
0.001210447933812 
-0.000962891786260 
0.000302228324383 
0.000000103004264 
0.000000000189553 
0.000000055144587 
0.000000112101364 
0.000000003575231 
-0.000000000398924 
0.000000001344731 
-0.000000000302993 
-0.000000001388141 
0.000000000000183 
0.000000000000440 
0.000000000000000 


0.054100841628481 

0.006089078688449 

0.008467361236503 

-0.021693066305532 

-0.000804346804545 

0.000176203862846 

0.000052121483427 

0.001252911553069 

0.000107558177780 

0.006296309320745 

-0.010064387555574 

-0.000001070541827 

0.007016262114685 

0.000000532027276 

-0.009464088570511 

0.000000000000000 

-0.169040041511134 

-1.928648189364952 

-5.541468744601931 

1.997318957894570 

0.872184481452392 

-0.141216029956213 

-0.747218737106791 

-0.033288989432478 

0.000046090125425 

0.000031936036514 

0.000078817300583 

-0.000094878787391 

-0.000000054649459 

0.000001718567633 

0.004219981300979 

0.275968586078390 

-0.101363240880644 

-0.001135783571085 

0.325185108704812 


-0.041785235726328 

-0.035662322704087 

0.015713341388969 

0.013516676525586 

0.000711298455703 

0.010710860308677 

0.003122205523126 

-0.000064404164098 

-0.000314763306762 

-0.000001303184149 

-0.000000139545383 

-0.006360025457770 

0.000998606118169 

-0.000151815242906 

-0.000000656704239 

0.000000000000000 

-0.000117535002353 

0.000599141638988 

0.001578138953716 

-0.000835019645022 

-0.000196352161334 

0.001354586899106 

0.000653421916173 

0.000001204490578 

0.000000041073598 

-0.000000013011312 

-0.000000028540982 

-0.000000083473916 

0.000000001688860 

-0.000000000085675 

0.000000013441555 

0.000000108444688 

0.000000000000800 

-0.000000000041803 

0.000000000000000 


-0.033880893331203 

0.082310459499198 

-0.065476641465701 

0.020551526058053 

0.000007004289984 

0.000000012889615 

0.000003749831911 

0.000007622892726 

0.000000243115696 

-0.000000027126812 

0.000000091441720 

-0.000000020603548 

-0.000000094393576 

0.000000000012456 

0.000000000029921 

0.000000000000000 

-0.000000000002061 

-0.000001344987463 

-0.000003784607242 

0.000001520058881 

-0.000000611900468 

0.000000280357456 

0.000000852496513 

0.000000002162050 

0.000000000006158 

-0.000000000033081 

-0.000000000062322 

-0.000000000015207 

-0.000000000000018 

-0.000000000000707 

0.000000000026640 

0.000000000195530 

0.000000000000912 

-0.000000000000369 

0.000000000000000 


7 - AN ELEMENTARY CONTROL STUDY 

The purpose of this report is to present the details of generating the 
mathematical model to be used in investigating the control problem of a CGRC. A 
commensurate presentation of the full investigation of the control problem will 
require a report of similar magnitude. However, it is desirable to present a 
sufficiently detailed control study to illustrate the use of the model. 
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STATEMENT OF THE CONTROL PROBLEM TO BE INVESTIGATED 

The positive eigenvalues of the matrix A reflect a clearly defined problem 
with the model behavior. If intervention via the control variables is not invoked, 
the eigenstates will tend to go to infinity with time, and thus all descriptive states 
which are related to the canonical states (via the mapping MR) will tend to go to 
infinity with time. Thus, it is necessary to manipulate the control variables in such 
a way that the resulting controlled system will have no eigenvalues with positive 
real parts. In the following, we demonstrate a procedure for choosing feedback 
control so that the canonical state having the most positive eigenvalue (eigenstate 
15). Left to itself, this eigenstate ( 915 (b) w dl tend to infinity according to e 
(coefficient of t is rounded). The Science Advisory Working Group did not address 
the question of dynamic properties requirements, so we have no guidance from that 
source to help with the choice of parametrization of the dynamic properties of the 
controlled system. Arbitrarily, we select the criterion that the dynamics of 
eigenstate 15 should decay according to e~ u and ask what feedback control v= Kz 
will provide the desired behavior. The first question is “does a K exist such that we 
can obtain this behavior?” If the answer to the first question is “yes,” then the 
second question is “will the gain matrix K have impractically large elements? (as a 
rough rule of thumb large gains are difficult to attain and are thus undesirable). 

The procedure we will demonstrate is due to C. Blackwell (1991) and 
distributes the burden of assigning the new eigenvalue among the control variables 
according to the control authority of each individual control variable. In this way, 
the wasteful assignment of load to a control variable with little authority is 
avoided, advantage of a control variable with high authority will be taken, and 
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those control variables having equal authority will be weighted equally. The 
fundamentals are relatively straightforward. 

For the case when all disturbances are zero, the dynamic equation of 
canonical state 15 is 


015 = A 15 015 + ♦'jS* 


where ir 15 is the fifteenth row of BR. By using a control law which relates « to j J5 , 

u = *15 015 

we make the dynamics of ? 15 dependent on ? 15 only 


or 


015 - ^15 015 + * r 15*l5 015 


015 ~ (-^15 + * r 15*is) 015 

We now select k 15 to be related to br 15 in a very useful way: 

*15 = ° * r i5'/ll* r i 5 || 2 

and this choice gives us a very simple expression for the dynamics of ? 15 

?15 = (^15 + a ) 0|5 

Now suppose we would like the controlled system to exhibit dynamics 


015 ~ ^ 15 ?15 


then we must have 


which tells us that 


(Aj 5 + o) — A 1S 


o — A is - Ajj 

Now we know ^ 5 . In terms of the entire canonical state, what we have done is 
represented by 


0 = (ER+BRK q )i+ CRd 

where 
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k, = [0|*, 5 |0] 

(k 15 is in the fifteenth column). K q is the canonical state gain matrix. To find the 
gains on the descriptive states, one uses the fact that q = MR~ l x. Thus 


K — K q MR~ 1 


NUMERICAL DETAILS 

The following is an edited MATLAB dialog transcript of the process of 
finding K. The sequence of MATLAB transactions follows the sequence which was 
followed in the preceding explanation. The object is to change an eigenvalue so 
that it has a desirable value. 


Iml5=er(15,15)//It is also equal to DBER(1,15). 

LM15 = 

0.024455993417775 

<> 

Li5 = -l.0//We want to replace LM15 with this value. 

L15 = 

-1 

<> 

alph=115-lml5 
ALPH = 

-1.024455993417775 

BR15=BR(15,l:ll);K15=ALPH*BR15V(NORM(BR15)**2) 

K15 = 

0.004845071804937 

0.000998659667081 

-0.000024779501890 

-0.001035565963037 

-0.001052970165541 

-0.001219941869759 

1.569397322971582 

-0.016593198483799 

0.000091585233884 

0.004694007912200 

0.000209987492911 

<> 

KQ(11,19)=0;KQ(1:11,15)=K15 
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KQ = 

Columns 1 thru 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 5 thru 8 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 9 thru 12 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 13 thru 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 17 thru 
0.000000000000000 
0.000000000000000 


0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 

0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 

0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 

16 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

19 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.004845071804937 

0.000998659667081 

- 0.000024779501890 

- 0.001035565963037 

- 0.001052970165541 

- 0.001219941869759 

1.569397322971582 

- 0.016593198483799 

0.000091585233884 

0.004694007912200 

0.000209987492911 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 
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KQ(continued) 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

<> 

=KQ*MRNV 

Columns 1 thru 1 
0.000000001761492 
0.000000000363076 
-0.000000000009009 
-0.000000000376494 
-0.000000000382822 
-0.000000000443526 
0.000000570575775 
-0.000000006032683 
0.000000000033297 
0.000000001706571 
0.000000000076344 

Columns 5 thru 
-0.003081479818114 
-0.000635150464878 
0.000015759835571 
0.000658622976893 
0.000669692100514 
0.000775886592043 
-0.998141280880727 
0.010553322696618 
-0.000058248476228 
-0.002985402741147 
-0.000133552658766 

Columns 9 thru 1 
0.000001569924753 
0.000000323590773 
-0.000000008029180 
-0.000000335549338 
-0.000000341188736 
-0.000000395291755 
0.000508524084486 
-0.000005376612375 
0.000000029675912 
0.000001520976263 
0.000000068041213 

Columns 13 thru 
0.000000779273095 
0.000000160622719 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000502248067 

0.000000103522694 

-0.000000002568684 

-0.000000107348461 

-0.000000109152609 

-0.000000126461169 

0.000162686293142 

-0.000001720078092 

0.000000009493875 

0.000000486588537 

0.000000021767647 

8 

-0.000000004235251 

-0.000000000872964 

0.000000000021661 

0.000000000905225 

0.000000000920439 

0.000000001066395 

-0.000001371866426 

0.000000014504709 

-0.000000000080058 

-0.000000004103200 

-0.000000000183558 

12 

-0.003141769809613 

-0.000647577356628 

0.000016068181044 

0.000671509114729 

0.000682794808768 

0.000791067024432 

-1.017670186760809 

0.010759801328062 

-0.000059388123522 

-0.003043812958480 

-0.000136165652891 

16 

0.000000800529824 

0.000000165004128 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000523727462 

0.000000107949998 

-0.000000002678537 

-0.000000111939380 

-0.000000113820685 

-0.000000131869472 

0.000169643817340 

-0.000001793639821 

0.000000009899895 

0.000000507398228 

0.000000022698573 

0.000001551717042 

0.000000319837824 

-0.000000007936059 

-0.000000331657696 

-0.000000337231689 

-0.000000390707231 

0.000502626311729 

-0.000005314255372 

0.000000029331736 

0.000001503336249 

0.000000067252083 

0.000000063887918 

0.000000013168492 

-0.000000000326747 

-0.000000013655144 

-0.000000013884639 

-0.000000016086355 

0.000020694332727 

-0.000000218800660 

0.000000001207658 

0.000000061895965 

0.000000002768930 

0.000000797670602 

0.000000164414789 


0.000000520498212 

0.000000107284390 

-0.000000002662022 

-0.000000111249173 

-0.000000113118879 

-0.000000131056378 

0.000168597811167 

-0.000001782580424 

0.000000009838853 

0.000000504269663 

0.000000022558616 

0.000001573623857 

0.000000324353227 

-0.000000008048099 

-0.000000336339970 

-0.000000341992656 

-0.000000396223154 

0.000509722284378 

-0.000005389280912 

0.000000029745836 

0.000001524560034 

0.000000068201534 

0.000000000729043 

0.000000000150269 

-0.000000000003729 

-0.000000000155823 

-0.000000000158442 

-0.000000000183566 

0.000000236148934 

-0.000000002496797 

0.000000000013781 

0.000000000706313 

0.000000000031597 

-0.003163812931494 

-0.000652120855186 
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K(continued) 

-0.000000003985493 -0.000000004094208 -0.000000004079585 0.000016180917780 
-0.000000166558665 -0.000000171101992 -0.000000170490874 0.000676220522043 
-0.000000169357928 -0.000000173977612 -0.000000173356223 0.000687585398181 
-0.000000196213372 -0.000000201565609 -0.000000200845685 0.000796617267732 
0.000252419191705 0.000259304590911 0.000258378442745 -1.024810311378582 
-0.000002668821775 -0.000002741620928 -0.000002731828787 0.010835293686337 
0.000000014730413 0.000000015132224 0.000000015078177 -0.000059804799384 
0.000000754976236 0.000000775570204 0.000000772800130 -0.003065168800599 
0.000000033774031 0.000000034695306 0.000000034571386 -0.000137121011261 
Columns 17 thru 19 

0.000000182112913 0.000000289450152 0.002044610915727 
0.000000037536868 0.000000059661075 0.000421432444888 
-0.000000000931393 -0.000000001480356 -0.000010456901794 
-0.000000038924074 -0.000000061865899 -0.000437006830286 
-0.000000039578250 -0.000000062905646 -0.000444351370026 
-0.000000045854257 -0.000000072880728 -0.000514813105746 
0.000058989325709 0.000093757597701 0.662282630030616 
-0.000000623692659 -0.000000991296726 -0.007002297615535 
0.000000003442437 0.000000005471407 0.000038648791279 
0.000000176434837 0.000000280425421 0.001980862245636 
0.000000007892852 0.000000012544894 0.000088614315217 
<> 

//Check to see that the K just computed does generate the desired eigenvalue 

AC=A-fB*K;EIG(AC)//AC is the state coefficient matrix of the controlled system 
ANS = 

1.0D+04 * 

-5.266821895634869 O.OOOOOOOOOOOOOOOi 

-0.021360198075745 O.OOOOOOOOOOOOOOOi 

-0.0039772634929 1 8 0.009552095834876i 

-0.003977263492918 - 0.009552095834876i 

-0.000459990795721 - 0.000120566056908i 

-0.00045999079572 1 0.000 1 20566056908i 

-0.000021248217594 O.OOOOOOOOOOOOOOOi 

-0.000093892720314 O.OOOOOOOOOOOOOOOi 

-0.000100000000000 - O.OOOOOOOOOOOOOOOi 

-0.000366233205798 0.000176262166536i 

-0.000365106096282 0.0001 737867 17458i 

-0.000365101237584 0.000173829082129i 

-0.000366233205798 - 0.000176262166536i 

-0.000365106096282 - 0.000173786717458i 

-0.000365101237584 - 0.0001 73829082129i 

0.000000294968295 O.OOOOOOOOOOOOOOOi 

-0.000000000002605 O.OOOOOOOOOOOOOOOi 

0.000000027480838 O.OOOOOOOOOOOOOOOi 

0.000000000000000 O.OOOOOOOOOOOOOOOi 

<> 

//The result confirms the correct execution of the procedure. 

<> 

exit 

(END OF MATLAB SESSION) 
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One can easily confirm that the eigenvalue .000002445599342 does not 
appear in the list of eigenvalues of the controlled system, and that an eigenvalue 
with value -1.0, the value which was desired, has appeared. All the other 
eigenvalues are unchanged. No element of the gain matrix is appreciably greater 
than 1.0 which implies, as a course rule of thumb, that one may expect to 
encounter no difficulty in accomplishing the eigenvalue change. The basic 
procedure is seen to be simple, but a negative side is that to change other 
eigenvalues by this method will require successively obtaining the RJCT after each 
reassignment and repeating the process we have demonstrated. 

Several matters of interest are spawned by this demonstration. One is that 
literally all 19 states must be known and 165 gain elements are required to set this 
eigenvalue. However, good results can often be obtained by setting relatively small 
gains to zero, thus not requiring full state feedback. In the present example, it was 
found that setting all gains to zero except those having order of magnitude one 
reassigned the target eigenvalue to a stable value, but it was not -1.0. Also, many 
other eigenvalues were changed, but none changed to unstable values (compare to 
root locus principles). Another matter of interest is the behavior of the volume 
change descriptive state (the nineteenth state). This state is not controllable since 
its only influence is the glove manipulation disturbance (note that this state is 
connected only to the nineteenth canonical state and that the last row of BR is 
zero). This might appear to be a problem for system control, but it is not in fact. 
If, for example we replace V with P as a state, we find that all nineteen states are 

controllable. 
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8 - CONCLUDING REMARKS 


The principal thrust of this study was to begin familiarization with the 
control problem specific to a CGRC system. It was brought to the point of 
displaying the system characteristics which are relevant to the control problem. A 
number of processes to make the control problem study meaningful are illustrated, 
but not all the actions which should be taken were done in this study. An 
important example is that the system linearized state variable format model 
matrices should be balanced to maximize numerical effectiveness and to make 
physical judgment as effective as possible. The modal matrix is normalized, 
however. A serious study should include both these important enhancements in 
preparation for the control law synthesis. 

The control synthesis demonstration is very elementary, but it suffices to 
demonstrate that the local stability problem is resolvable via straightforward 
multivariate means. Although canonical space eigenvalue placement was used in 
the demonstration, alternate choices such as linear quadratic regulator synthesis or 
minimum norm gain matrix synthesis could have been used. To describe the full 
control study will require an additional report which is based on the contents of 
this report but which is dedicated to the synthesis activity. The uncertainties in 
the model are yet to be characterized, and as a consequence, robustness 
considerations cannot be brought into the control design process at this stage. 
Because of the dynamic dimensionality and the hypothetical nature of the object of 
the study , it is not certain that any significant constructive purpose can be served 
by further addressing the control problem via this model. Stated differently, the 
resource might be better expended in the design of the control system for a model 
which is generated from a more completely defined physical system. 
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ABSTRACT 

The Crop Growth Research Chamber (CGRC), which 
is being developed unoer the Controlled Ecological Life 
Support System program at NASA Ames Research 
Center, will be operated to support research on the growth 
dynamics of crops of higher plants within a cloned, 
precisely controlled environment. The CGRC is the first 
in a senes of instruments which will be incorporated into 
bioregenerative life support systems for space habitats. 
The dynamic processes of the thermal and fluid portions of 
the CGRC are profoundly coupled with thoee of the plants 
and there is strong reason to believe that consideration of 
these interactions is necessary in order to design a CGRC 
which will support sound research. 

In this paper we aescnbe the modeling system which 
we have deveiopeo for the thermal and fluid dynamics of 
the CGRC. The bases of the modeling system are 
symbolic representations of the individual processes which 
are included in the representation. The system includes 
the plant growth chamoer. the associated control system 
components and the control devices. An example of the 
derivation of the dynamic equations of the system from 
the symbolic modeling system is presented. 

AN IMPORTANT ELEMENT of the Controlled 
Ecological Life Support System (CELSS) Program at 
NASA Ames Researcn Center is the development of a crop 
growth researcn chamoer (CGRC). This CGRC contains a 
plant growth chamber (PGC) requiring a precisely 
controlled internal environment. The CELSS program and 
• ts recent status are described in (l). 

The CGRC is one of a sequence of physical devices 
which are to t>e developed under the CELSS program, 
leading eventually to fully cycling bioregenerative life 
support systems lor space or planetary habitats i2|. It will 
be used to conduct researcn on the growth of crops of 
nigher mants from seed to msturation employing stringent 
control of all environmental parameters (e.g., lighting, 
temperature. atmosonertc oxygen. carbon dioxide 
concentrations, humidity, and pressure, nutnent solution 
temperature, ion ana dissolved oxygen concentrations, and 


pressure) which affect the production of biomass, oxygen, 
sad transpired water vapor, and the consumption of 
carbon dioxide, liquid water, and nutrient# by the plants. 
Knowledge of the au (ecology of the crop plant populations 
to be used in a CELSS will permit utilisation of the 
CELSS plant growth unit ns a multiple input-multiple 
output (MIMO) life support system component whose 
function and performance can be vaned as needed by 
application of the appropriate control strategy. 

The prototype chamber unit currently under 
develop meat is a precursor to a specs- based unit (CROP) 
that will be capable of conducting experiment# on the 
* fleet of the space environment on the growth of plant 
crape. It is desirable to regulate environmental conditions 
within the space-based unit to the same degree of accuracy 
as those within the ground-based unit in order to separate 
the effects of flight mission variables, principally gravity, 
from the effects of controllable environmental parameters 
on plant growth response. Suggested ancillary usee of the 
chambers are to “develop' control systems, to validate 
models of plant growth and system operation, to explore 
the development of expert systems for operation of plant 
growth devices, and to explore the use of automation ana 
robotics and other devices to save human labor” (3). 

Both the ground-b a se d and the space-based unit* will 
be designed to provide environmental conditions that are 
independent of thoee of the ambient environment. While 
the space- ba s ed unit and follow-on versions of the 
prototype ground-based unit will be closed with respect to 
the exchange of materials, the initial ground based unit 
will emulate closure by including sufficiently large storages 
and sinks of materials which require removal or 
replenishment during particular pnases of plant growth. 

Tolerance specifications for environmental vanables 
for the ground-based and space-based units were 
developed and modified at workshoos held in Seotember 
1984, April 198ft and June 1988 (4. 5). The tolerances on 
environmental variables are described in i6) and are given 
in Table 1. The CGRC differs from other plant growth 
chambers described in the literature i7 * 12) in size, degree 
of material closure, tightness of control, and number of 
controlled environmental vanaoies. 
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TABLE 1. Tolerance requirements for the 

atmospheric variables in 

the shoot tone of tbe plant 

growth chamber of tbe 
Research Chamber. 

NASA CELSS Crop Growth 

Air temperature 

8-40* C ± 1* C 

Relative humidity 

38-90% ± 2% 

Carbon dioxide 

23-3000 ppm ± 0J2% 

Oxygen 

8-25% ± 5% 

Nitrogen 

78-06% ± 5% 

Gage pressure 

(US* HfO ± 0.25" 

Air velocity 

0.5 m see* 1 ± (•) 

Photosynthetic 

0-3000 pmolee nT 1 s~ l 

photon flux 

± 10 pmoles s" 1 

Surface temperatures 

Air temperature + T C 
• not established 


The dynamic processes of the thermal and fluid 
upt cU of the CGRC an profoundly coupled with those of 
the pUaU and there Is strong roaaoo to bettave that 
c ooa ldar mtl on of Hum interactions Is nsessssry In ardor to 
4mixn a CGRC wkick will support sound rasuarck. For 
tko cans of complex systems, practical nsosasity prodkatos 
tko systematic examination of control Issues on tke use of 
adequately accurate mathematical models of tke system 


Many very Important questions related to the 
acceptability of performance of a passively controlled 
system with closure can be evaluated by analyrinf a 
mathematical model of tke dynamic behavior of the 
system. Due to disturbances, uncertainties, or variations 
in physical or- biological characteristics of the system 
c om ponents, passive control may foil to provide acceptable 
performance, la that case, recourse to the use of active 
control Is a compelling alternative since active control Is 
well established to have highly reduced sensitivity to 
uncertainty. The same mathematical models used to 
determine the adequacy of passive control can and should 
be need to study the requirements for active control laws 
which will provide stable, desirable system performance in 
which compensation is provided for disturbances and for 
uncertain or variable characteristics of elements of the 


oontrouea system. 

Uncertainties in the CGRC system are most manifest 
In the living component. The inherent vnrinhUity of 
biological systems (13) makes the teak of proems 
mathematical description of their behavior difficult. 

The term “stable" is used hero in the sense of 
stability to a region, i-e. the states of the system remain 
within an acceptable tolerance of the desired system state. 
The region of tolerance is described by the specifications 
on the environmental variables given in Table 1. 

For a plant growth unit operating as a regulating 
system, this implies prescribed bounds on the deviation of 


the actual system variables from dsewad 
a plant growth unit operating aa the M1MO tracking 
system described earlier. U., one which is tranaitioeieg 
bom an initial to a final state along a prescribed 
performance path, this implies proKribed boned. on the 
deviation of tke actual system variables from the <M»rod 
trajectory between initial and final states- A control 
which provides stability aa defined is said to be a 
control for . the system subject to the specified 


Control strategies, proposed or implemented, hove 
men reported for growth chum bars, greenhouses, end 
ivestock bsildiag. |«, 11. 12, If - 

itreteciss heve been developed using ed hoc methods or, 
to varying degress. etUiring control theory formalom. 
Most depend apoa some farm of episodic manual toeing to 
retain darirobla performance. It is possible that * 
control may achieved by manual tuning, bat it is 
fcA know whether robustness bus been achieved or wbetner 
future aystem partnrbnthma will require “ 

desirable to derrinp a control syatem that requires little or 
no human intervention, ea will be the case for the apace 
home system, a formal procedure for the synthesis of a 
control strategy is required. ■ 

Thu robustness objectives are moat readily 
through the development and utilisation of mathematic^ 

■ of the r of the CGRC. i-e- the PGC and 

mats of the control syatem. The 
[ which will be develo p ed for the 
;GRC is baaed upon the philosophy that a modal which to 
toea topad for control system design to one in which 
Jophyrical and physical phenomena am mpraasatad in the 
nost pareimoaious format which to adequate far the 
analytical and computational Beads of control design is 
oa treat to the repre se ntation of mterooeopic »“• f " 
■sue of control modeling and scientific research modaiiag 
s discussed la {3}. 

Approaches to robust control theory under current 
levelopment full In five major categories which cun be 
Minfiy described am (1) H« methods [24], ( 2) Baser 
jptimai quadratic Gnauaiaa with loop treaefor recovery 
[LQG/LTR) [251, (2) quantitative ^backtheo^ (QfTT) 
“18], (4) variable structure with sliding mode. (VS/SM) 
17], and (5) Lyapunov theory baaed methods [2*]. The 
lutbemuticul modeling structure developed in this P*P*r 


uuucui mooeung — T . , 

LB omsted toward the atilisutiou of (5), a peth °° 0 * 0 ^. 
which appears to be most promising for the types of 
mteme. uncertainties, end performance measures 
^counterod in CGRC research and CELSS application.. 
Illustrations of C ELSS minted applications of Lyapunov 
theory kneed methods of robust control design can be 
found in (20) and (30). 

The purpose of this paper is to illustrate the approach 
to the derivation of a continuous time state variable 
format model for use in developing a robust 

control law for tbe CGRC. 

SYSTEM DESCRIPTION 

The development of a system to «ont*ol the 
environmental conditions within the PGC of the CGRC 
involves (1) the selection of component devices (*•*• 
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actuators) which will affect the environmental variables 
that are to be controlled, (2) the selection of component 
devices which will provide suitable measures of the 
environmental variables (i.e. measurement instruments), 
and (3) the development of control algorithms which 
utilise the output# of the measurement instruments to 
provide the necessary combinations and time sequencing of 
inputs to the actuators so that the required pe r f or mance 
characteristics of the environmental variables within the 
CGRC are met. 

Fig. 1 is a representation of the arrangement of the 
components hypothesized for control of the atmosphere of 
the plant growth chamber. The plants receive radiant 
energy from a light source above the chamber (not 
shown). Water and nutrients are supplied to the roots by 
means of a nutrient delivery system that is materially 
isolated from the chamber atmosphere. Air flow into the 
chamber is assumed to be sufficient to assure that uniform 
conditions exist in the atmospheric control volume within 


the chamber and surrounding the plant canopy control 
volume. Gaseous exchange of carbon dioxide, water vapor 
and oxy gen occurs betwee n the canopy and the 
atmosphere. 

Air enters near the top of the upper chamber and is 
thoroughly mixed with the air in the upper portion of the 
chamber. The resulting mixture flows b e twe en the walls 
and the baffle formed by the plant support surface into 
the lower portion of the chamber and then into the 
ducting located near the bottom. A HEPA filter is 
provided to remove particulates from the air as it leaves 
the chamber. Air flow out of the filter is affected by the 
controllable orifice flow area of a valve. A portion of the 
air flow is diverted, by means of a controllable flapper 
valve and fan, into a gas separator which removes excess 
oxygen or carbon dioxide. A centrifugal pump (blower) 
serves to compensate for pressure losses within the system 
and provide the required air movement within the system. 
Makeup gases are injected into the flow it ream to 
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maintain the required atmospheric composition. A portion 
of the flow is diverted through a dehtanudifying heat 
exchanger, where the condensate is removed from the 
system. Two variable flow area orifices are present, one 
each in the flow path through the dehumidifier and in the 
parallel bypass. The orifice flow areas are variable in order 
to regulate the mass flow ratio of the paths. The air in 
the two flow paths is mixed and flows through a section of 
ducting. A portion of the flow is diverted through either a 
beater or a cooling heat exchanger. Three variable flow 
area orifices are present one each in the flow path through 
the cooling heat exchanger, the heater and the parallel 
bypass. The orifice flow areas are variable in order to 
regulate the mass flow ratio of the paths. Since meeting 
performance specification is considered paramount to 
economic constraints in this design concept, two 
independently controlled heat exchangers are utilised in 
humidity control and temperature regulation. The flows 
are mixed and flow through the ducting to the chamber 
inlet. 

The task of modeling the processes which affect the 
dynamics of the thermal environment of greenhouses, 
plant growth chambers, livestock houses and confuted 
spaces has been addressed (9, 11, 12, 14, 18, 19, 31 - 37]. 
The processes included in the model depend upon the 
assumptions of importance of each process in determining 
the dynamics of interest for the particular system being 
modeled. While the physics of the processes which have 
been included in these models are applicable to the 
dynamics of the CGRC, they do not indude all the 
p r o cs ss ss necessary to account for the dynamics of the 
variables of interest in the CGRC. 

In the following sections, we describe the assumptions 
which we have mode concerning the function of each 
component of the CGRC and the physical and biophysical 
processes which we hypothesise to contribute significantly 
to the thermal and fluid dynamics within them. In this 
formulation, we consider the scenario in which the PGC 
contains a crop of plants forming a dosed canopy and the 
system goal is to maintain temperature, relative humidity, 
pressure, carbon dioxide concentration, oxygen 
concentration, and mean air velocity in the chamber 
within set tolerances about constant operating points. 

MODELING PROCEDURE 

The process of mathematical model development is en- 
hanced by the preliminary formulation of a s y mb olic 
model of the system. The symbolic model is a schematic 
rep res e n tation of each of the generic elemental physical 
and biological processes which are to be ioduded in the 
mathematical model and the manner in which these 
processes interact in the particular system being analysed. 
The mathematical equations describing these elemental 
processes and the structure of their interactions combine 
to become the mathematical model of the dynamic 
behavior of the system. Other symbologies have been used 
for developing models of engineering systems for control 
system studies. These indude linear graphs and dreuit 
diagrams (38, 39] and bond graphs [40, 41). In this paper 
we present the symbology of the modding system which 
we have developed for the thermal and fluid dynamics of 


the atmospheric control system of the CGRC. 

ASSUMPTIONS — Some sssumptions are made in 
order to achieve the minimum complexity model for the 
current analysis. These assumptions may be modified u 
more spedfic details of the system become available. 

(1) Leakage from the system is negligible. 

(2) The air in the system consists solely of nitrogen, 
oxygen, water vapor and carbon dioxide and behaves as 
an ideal gsa. 

(3) The boundary of the system indudes the PGC and 
the elements of the environmental control system. 

(4) The air ducts of the system are perfectly insulated. 
Thermal capacitance of the air duct walls is negligible. 

(5) Distributed phenomena are adequatdy represented 
by lumped models. 

(6) The radiant energy source consists of lamps (and 
their associated reflectors) which produce 
photosyntheticaiiy active radiation (PAR) as well as long 
wave radiant energy. Control of the radiant energy 
intensity at the canopy level, frequency content, and 
sequendng (Le. effective photoperiod) is unaffected by 
PGC processes. A substantial portion of the long wave 
radiant energy produced by the radiant energy source in 
the 2800 nm to 100,000 nm range is absorbed by a liquid 
water/giaas filter [42]. The water in the filter is cooled by 
a separate control system which maintains the water at a 
prescribed temperature. 

(7) The control system for the nitrogen gas blanket 
which separates the inner walls from the exterior walls of 
the PGC is able to maintain the temperature of the inner 
walls at a prescribed temperature above the dewpoint of 
the chamber air. 

(8) The canopy processes of transpiration, photo- 
synthesis, photorespiration and maintenance respiration 
are the only canopy processes of significance for the time 
scale considered in this model. (See Fig. 3.3 in (43).) 
Growth is negligible, i.e. the leaf carbohydrate pool as 
defined in (44] is an infinite source/sink. Similarly, there is 
an infinite water source for internal water associated with 
the plant. 

(9) The plant population can be modeled by ag- 
gregating the sum of individual plant processes and their 
interactions into a single canopy process. The canopy 
biomass is sufficiently large to affect the fluid and thermal 
dynamics of chamber. The conditions within the chamber 
affect the physiological dynamics of the canopy. 

(10) A nutrient delivery system in the lower chamber 
is maintained at fixed temperature above the dew point of 
the chamber air. There is no material interchange 
between the nutrient delivery system and the chamber air. 
Water and' nutrients are delivered to the plants at the 
ideal concentration, pH and pressure. Dynamics of 
nutrient delivery and plant uptake and translocation are 
negligible. 

(11) The flow rates and entrance temperatures of the 
working fluid in the beat exchangers are maintained by 
separate control systems whose dynamics are negligible. 

(12) An infinite sink exists for the removal of gases 
and condensate. 

(13) An infinite source exists for the addition of gases 
and water vapor. 

(14) The dynamics of air diversion from the main flow 
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Fig. 2. Conceptual model of the Plant Growth Chamber and atmospheric control system of the Crop j 
Growth Researcn Chamber indicating points and flow path designations for symbolic modeling. Points are I 
indicated as encircled numbers. Flow paths are shown as numbered arrows. Arrows indicate the direction 
of assigned positive flow. 


path, gas separation, component gas removal, and reintro- 
duction of the diverted gas into the mam flow path can be 
neglected. The process can be represented by 
instantaneous removal of the gas component from the 
main flow stream. 

(15) The dynamics of source gas or water vapor in- 
jection are negligible. The process can be represented by 
instantaneous addition of the gas or vtpor to the main 
flow stream at a prescribed temperature. 

(16) Flows are turbulent but fully developed. 

(17) Other than in the blower, incompressible flow 
relationships are adequate to describe system processes. 

ASSIGNMENT OF POINTS AND FLOW PATHS - 
Fig. 2 re present i the conceptual model of the CGRC with 
the points and flow paths needed for the modeling 
symbology assigned to the system. A point is the location 
at which ail the v&nabies representing potentials with 
respect to a fixed reference within a designated volume 
may be assumed to be concentrated. Point numbers are 


shown encircled. Flow paths represent the constitutive 
relationships (processes) which occur between points. 
Flow paths are shown which arrows indicating the flow 
direction which is assigned a positive value. Path 
numbers are not circled. 

CONSTITUTIVE RELATIONSHIPS - The symbol- 
ogy for the constitutive relationships is representative of 
processes between and at points. In many cases, the 
symbol represents aggregated phenomena, i.e. total molar 
flow rates. 

Storages — We account for molar storages and energy stor- 
ages in the system. Fig. 3 illustrates the molar storage at 
point 1. The symbol illustrates the rate of change of upper 
chamber atmospheric mass as a function of molar flow 
rate into the chamber through flow path 1 (from the 
atmospheric control system) (n^) and flow pmh 2 (from 
the plant canopy, through gaseous exchange of the 
products of transpiration, photosynthesis, photores pi ration 
and maintenance respiration) (n^) and molar flow rate out 
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of the upper chamber through flow path 3 (to the 
atmospheric control system) (n 3 ) end flow path 2 (to the 
plnnt canopy* through gaseous exchange of the reactants 
for photosynthesis* photorespiration and maintenance 
respiration) (n 2 ). Significant molar storages in the system 
are located at points 1 (upper chamber atmospher e)* 2 
(plant canopy), 5 (lower chamber atmosphere), 15 (blower 
entrance), and 36 (gas resupply chamber). Within the 
plant canopy, molar storage is composed of three 
substorages: substomatal cavity (gaaeous products) (point 
2), leaf carbohydrate pool (point 46), and internal plaat 
water pool (point 44). (Points 46 and 44 are not shown 
on Fig. 2.) The remaining small amounts of molar storage 
in the system are lumped into these storages. 


Fig. 4 illustrates the energy Stonge at point 1. The 
symbol illustrates the rate of change of upper chamber 
atmosphere energy ns a function of enthalpy flow rate into 
the upper chamber atmosphere through flow path 1 (from 
the atmospheric control system) (&t) end flow path 2 
(from the plaat canopy, through gaaeous exchange of the 
products of transpiration, photosynthesis, photorespiration 
and maintenance respiration) (&?), enthalpy flow rate out 
of the upper chamber atmosphere through flow path 3 (to 
the atmospheric control system) (H§) and flow path 2 (to 
the plant canopy, through gaaeous exchange of the 
reactants for photosynthesis, photorespiration and 
maintenance respiration) (H$), hsat transfer into the 
upper chamber atmosphere from the plant canopy (Qj) 
and from the walls and lights ($41), *®d work (PjVj) 
done at the moving boundary (glove ports). These 
storages are located at points 1 (upper chamber 
atmosphere), 2 (plant canopy), 5 (lower chamber 
atmosphere). 15 (blower entrance), and 36 (gas resupply 
chamber). An additional energy storage within the plant 
canopy represents energy storage associated with the 


chemical bonds in the leaf carbohydrate (point 46). 
Mechanical energy storage occurs at point 37, is illustrated 
in Fig. 14. and is discussed in conjunction with blower 



Inortnnee - We account for total fluid inertia within the 
system by lumping all inertance effects into two fluid 
inert sn ret one in each section of the control system 
between the blower and the PGC. Fig. 5 illustrates fluid 
inertance located between points 13 and 14. The rate of 
change of the molar flow rate between points 13 and 14 
( n ) is proportional to the total pressure difference 
between points 13 and 14. A similar process is located 
between points 35 and 3. 



* *1 + *•* 


~ Dj 


Fig. 3. Symbology and sample equation for molar 
storage processes The process is illustrated for the 
molar storage at point 1 on Fig. 2. Molar flow rate n< 
represents the sum of all atmospheric species in that 
flow. Flow path arrows indicate the direction of 
assigned positive flow. 
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Fig. 5. Symbology and sample equation for fluid 
inertance processes. The process is illustrated for fluid 
inertance between points 13 and 14 on Fig. 2. Arrows 
indicate the direction of assigned positive flow. 
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Pipe Friction - We account for total pressure losses due 
to friction in the ductwork (constant area) within the 
system by lumping all constant area friction effects into 
two friction loss processes, one in each section of the 
control system between the blower and the PGC. Fig. 6 
illustrates pipe friction located between points 12 and 13. 
The total pressure difference between points 13 and 14 is 
proportional to the square of the molar flow rate between 
points 13 and 14 (n 8 ). The friction factor is determined 
by the duct surface and the Reynolds number. Effects of 
fittings and bends are included in equivalent length. A 
similar process is located between points 24 and 25. 



P 12 “ P 13 
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Fig. 6. Symbology and sample equation for pressure 
changes due to pipe friction for constant area, single 
phase flow. The process is illustrated for friction 
between point 12 and 13 on Fig. 2. Arrows indicate 
the direction of assigned positive flow. 


Flow Splitting and Joining - We account for total 
pressure changes due to flow division and flow merging. 
These processes occur in the section of the control system 
involving dehumidification and temperature balancing. 
Fig. 7 illustrates flow division between points 17 and 
points 18 and 19. The total pressure difference between 
point 17 and point 18 (or between point 17 and point 19) 
is a bilinear function of the molar flow rates (n 14 ) and 
(n l6 ) (or (n l4 ) and (n l5 )). The weighting of each term 
(a x 2 3 ) is a function of the splitter geometry. A similar 
process is located between points 25 and points 26, 27, 
and 28. 

Fig. 8 illustrates flow merging among from points 22 
and 23 to point 24. The total pressure difference between 
point 22 and 24 (or 23 and 24) is a bilinear function of the 
molar flow rates (n 21 ), (n 22 ), and (n 23 ). A corresponding 
process relates the flows through points 32, 33 , 34, and 
35. 

Sudden Expansion - We account for total pressure 
changes due to sudden increases in cross sectional flow 
area. Fig. 9 illustrates sudden expansion from point 3 to 
point 1. Pressure change varies as a function of the 
square of the molar flow rate between points 3 and 1 (n^. 
Similar processes are located between points 8 and 9, 14 
and 13. and 16 and 36. 

Sudden Contraction - We account for total pressure 
changes due to sudden decreases in cross sectional flow 
area. Fig. 10 illustrates sudden contraction from point 5 


to point 6. Pressure change varies as a function of the 
square of the molar flow rate between points 5 and 6 (n 5 ) 
and as a function of the upstream and downstream duct 
diameters. Similar processes are located between points 
10 and 11. and 13 and 17. 



n i4 = n i» + "is 


+ 


Fig. 7. Symbology and sample equations for flow 
splitting process. The process is illustrated for flow 
splitting between point 17 and points 18 and 19 on 
Fig. 2. Pressure relationship from [45]. Arrows 

indicate the direction of assigned positive flow. 


Flow through an Orifice — We account for total pressure 
changes due to flows through orifices. Fig. 11 illustrates 
flow through an orifice between points 11 and 12. 
Pressure drop across the orifice varies as a function of the 
square of molar flow rate between points 11 and 12 (n 7 ). 
For a circular variable diameter thin plate orifice Kj can 
be approximated as a quadratic function of the orifice to 
duct diameter. Similar processes are located between 
points 21 and 22, 20 and 23. 31 and 32, 30 and 33, and 29 
and 34. All of these orifice areas are shown as variable. 
These orifice areas are available as system control inputs. 
A similar process occurs between points 1 and 5 around 
the plant growth platform which forms a baffle between 
the upper and lower regions of the PGC. In that case, the 
orifice area is fixed. 

Flow through g, Porous Media — We account for total 
pressure changes through the particulate filter. Fig. 12 
illustrates flow through a filter between points 9 and 10. 
The pressure difference between points 9 and 10 is 
proportional to the square of the molar flow rate between 
points 9 and 10 (n 4 ). Filter porosity may decrease as 
material accumulates within the filter. 
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Fig. 8. Symbology and sample equations for flow 
merging process. The process is illustrated for flow 
merging from points 22 and 23 to point 24 on Fig. 2. 
Pressure relationship from [45). Arrows indicate the 
direction of assigned positive flow. 



Gas Removal — We account for the removal of excess 
atmospheric gases from the system. Fig. 13 illustrates 
selective removal of gases from the system at point 12. 
The molar and entropy flow rates downstream of the 
extraction point are equal to the upstream rates less the 
rate of gas removal. 
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Fig. 10. Symbology and sample equation for sudden 
contraction process. The process is illustrated for 
sudden cross sectional flow area contraction from point 
5 to point 6 on Fig. 2. The formulation of K, is 
specific to the sudden contraction flow process and 
may be approximated as proportional to the ratio of 
downstream to upstream duct diameters. Arrows 
indicate the direction of assigned positive flow. 
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Fig. 11. Symbology and sample equation for flow 
through an orifice. The process is illustrated for flow 
through the orifice from point 11 to point 12 on Fig. 2. 
The formulation of K< is specific to the orifice flow 
process and may be approximated as a quadratic 
function of the duct to orifice diameter. Arrows 
indicate the direction of assigned positive flow. 


QmA Injection — We account for the addition of 
atmospheric gases. This procsss is illustrated in Figa. 3 
and 4. In that case* the gases which are added to molar 
storage originate from another system component (the 
plant canopy). A molar flow rate of additional gases (n x7 ) 
from an external source is indicated at point 13. Molar 
and energy influx rates from the external sources are 
modeled in a fashion similar to that shown in Figa. 3 and 
4. respectively. 

Blower — We account for thermal, mechanical and fluid 
processes associated with the blower. The processes 
involved with the blower are illustrated in Fig. 14. These 
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processes occur between points 15. 37 end 16. Motor 
torque is availed ie as a control input. A mechanical 
energy storage (due to the rotational inertia of the motor 
and impeller) occurs at point 37. Mechanical power is 
added to the system by the pump motor. Pump efficiency 
is accounted for in the mechanical to fluid power 
conversion process equation. An isentropic compression 
process is modeled in the blower. 
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Fig. 12. Symbology and sample equation for flow 
through a porous media. The process is illustrated for 
flow through the porous filter from point 9 to point 10 
on Fig. 2. The formulation of K, is specific to the 
porous flow process. Arrows indicate the direction of 
assigned positive flow. 
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Fig. 13. Symbology and sample equations for gas 
removal. The process is illustrated for gas removal at 
point 12 on Fig. 2. Arrows indicate the direction of 
assigned positive flow. 


Heat Exchangers — The three heat exchangers (dehumid- 
ifying, cooling, and heating) are represented by a generic 
counterflow model. The processes involved in the heat 
exchanger model are illustrated in Fig. 15 for the 
dehumidifier between points 19 and 21 and 18. The 


processes in the working fluid, whose dynamics are 
separately controlled, occur between points 38 and 39. 
The logarithm mean temperature difference method is 
used. Since condensate is removed in the dehumidifier, a 
flow path (18) is shown for molar and enthalpy flow. The 
condensate flow path will not appear for the other two 
heat exchanger systems. For the dehumtdifier the molar 
flow rates for water vapor and other atmospheric 
component gases are considered separately. This process 
is not required for the other heat exchanger systems. 



Molecular Diffusion — We account for molecular diffusion 
of gases between the canopy and the chamber air. Fig. 16 
illustrates diffusion of water vapor between points 1 and 2 
(the substomataJ cavity). The diffusion resistance 
represents the sum of boundary layer and stomatal 
resistance. Boundary layer resistance various with wind 
speed [46]. The stomatal resistance is different for the 
different molecules in the atmosphere [47] and vanes as a 
function of canopy temperature, ambient relative 
humidity, canopy carbon dioxide concentration, and 
canopy water potential (48). Stomatal and boundary layer 
resistances are defined as overall canopy level resistances. 
The molar flow rate (nj) is proportional to the difference 
between specific humidity at point 2 and at point 1. 

We account for resistance to diffusion through the 
mesophyll cells for gases and. for carbon dioxide, the 
resistance associated with the carboxyiation reaction [47, 
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49, 50]. The process is illustrated in Fig. 17 for carbon 
dioxide. The molar flow rate (n„) is proportional to the 
difference between carbon dioxide concentration at point 2 
and point 40. 
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Fig. 15. Symbology and sample equations for the 
heat exchanger. The processes are illustrated for the 
heat exchanger between point 19 and 21 on Fig. 2. 
Arrows indicate the direction of assigned positive flow. 
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Fig. 16. Symbology and sample equation for molecu- 
lar diffusion of gases between the plant canopy and the 
chamber air. The process is illustrated for water 
vapor diffusion between point 1 and point 2 on Fig. 2. 
Arrows indicate the direction of assigned positive flow. 
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Fig. 17. Symbology and sample equation for diffusion 
through the mesophyll and resistance asso ci a t ed with 
the carboxylation reaction. The process is illustrated 
for carbon dioxide between point 2 and point 46. 
Arrows indicate the direction of assigned positive flow. 


Water Transport - We account for water transport 
through the plant from the point designated as an internal 
water pool (source) to water vapor in the substomatai 
cavity. The process is illustrated in Fig. 18. The molar 
flow rate (n^g) is proportional to the water potential 
difference between points 2 and 44. 
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Fig. 18. Symbology and sample equation for water 
transport through the plant. The process is illustrated 
for transport from point 44 to P 0 *® 1 2. Arrows 
indicate the direction of assigned positive flow. 


Convective H..t Tran.fcr - We account for convective 
heat transfer between surfaces in the chamber and the 
chamber air. Fig. 19 illustrates convective heat transfer 
between the plant canopy and the upper chamber air. 
The convective heat transfer rate (<i 2 ) i* proportional to 
the temperature difference between point 1 and point 2. 
A similar process occurs between points 1 and 4i and 
points 5 and 7i. 

Radiative Heat Trnnifcr ~ We account for radiative heat 
transfer between iutfaces in the chamber. rig. 
illustrates radiative heat transfer between the plant 
canopy and the upper chamber walls. The radiative heat 
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transfer rale (Q 40 ) »> proportional to the difference 
between the fourth power of temperature at point 2 and 
the fourth power of temperature at point 4i. 





UA 



Qi = be Ac (T 2 - Tj) 


Fig. 19. Symbology and sample equation for convect- 
ive heat transfer. The process is illustrated for 
convective heat transfer between the plant canopy, 
point 2, and the upper chamber air, point 1. Arrows 
indicate the direction of assigned positive flow. 


(without photorespiration) is illustrated in Fig. 22. 



f 4i — a A 1( K« 

Fig. 21. Symbology and sample equation for absorp- 
tion rate of photosyntheticaliy active radiation (PAR). 
Effective canopy area is used since ail leaves are not 
equally illuminated. Average energy content of photo- 
synthetic photons from [13, 42. 51). Incident PAR 
above the canopy includes PAR emitted from all 
sources as well as direct and indirect sources and 
scattered PAR (44). Arrows indicate the direction of 
assigned positive flow. 
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Fig. 20; Symbology and sampie equation for radiative 
beat transfer. The process is illustrated for radiative 
heat transfer between the plant canopy, point 2, and 
the upper chamber surfaces, point 4. Arrows indicate 
the direction of assigned positive flow. 


Absorption q! PAR - We account for absorption of photo- 
synthetically active radiation (PAR) by the plant canopy. 
Fig. 21 illustrates absorption by the plant canopy (point 
2) of PAR from the visible light energy at point 45. 

Bound Energy — We account for PAR energy absorbed 
but not used in photosynthesis, energy present in chemical 
bonds as a result of photosynthesis and the energy 
released as a result of respiration. The energy associated 
with the inorganic materials at point 2 are related to the 
energy associated with the organic materials at point 46 
through the photosynthetic and respiration processes. 
Some of the PAR absorbed is not used for photosynthesis 
and becomes part of the energy storage at point 2. The 
energy released as a function of photorespiration and 
maintenance respiration becomes part of the energy 
storage at point 2. Canopy data for wheat in a CELSS 
can be found in [13]. The process for gross photosynthesis 
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Fig. 22. Symbology and sample equations for rate of 
transformation of energy to/ from organic chemical 
bonds. The process is illustrated for photosynthesis. 
Average energy content of photosynthctic photons 
from [13, 42, 51). Incident PAR above the canopy 
includes PAR emitted from ail sources as well as direct 
and indirect sources and scattered PAR [44). 
Maximum rate and half saturation constants are 
functions of carbon dioxide concentration (44). Arrows 
indicate the direction of assigned positive flow. 


Modulated Sources — We account for modulated sources. 
Sources are devices which are capable of delivering energy 
to a system from a location external to it. Sources are 
considered to deliver a variable to the system in a 
prescribed manner. A dependent source is one in which 
that source variable is a function of another independent 
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variable. The second independent variable if termed the 
variable (38). No energy is drawn from the 
port of the modulating variable. A modulated source is 
illustrated in Fig. 23. An enthalpy flow rate. H 3 , occurs 
between point 1 and 2. The enthalpy flow rate if a 
function of the modulating variables n 3 , the molar flow 
rate in path 2, and k 3 , the specific enthalpy at point 2. 



Fig. 23. Symbology and sample equation for a modu- 
lated source. The process is illustrated for a modu- 
lated enthalpy flow rate source between points 1 and 
2. The enthalpy flow rate in path 2 is a flow source 
which depends upon molar flow rate in path 2 and the 
specific enthalpy at point 2. Arrows indicate the 
direction of assigned positive flow. 


STRUCTURAL RELATIONSHIPS - The tymbols 
for of the pro cents which were developed in the 

previous section are used in the formulation of the 
modeling diagram for the system. The modeling diagram 
describes tbe manner in which those processes which were 
previously described interact. Whereas the same process 
may occur at many points and is generic to many physical 
systems, the structure of the modeling diagram represents 
a proems interaction which is unique to the system being 
modeled. When the process symbols are United at the 
appropriate structural points, two modeUng diagrams are 
produced: one which represents the molar component and 
one which represents the energy component of the system. 

The molar component modeling diagram is shown in 
Fig. 24. The continuity principle is used to write the 
equations at point, i-a. the sum of the molar flow 

f£lai into a poiat is equal to the sum of the molar flow 
rates out of a point. Examples are: 

At point 35, 

“3 + ®31 + ®3J * *1 

At point 3, 

®i « n x 

and 

At point 1, 

Dl + n 2 = n 3 + hcnj. 

The energy component modeling diagram is shown in 
Fig. 25. The energy conservation principle it uted to write 
the equations at each point, i.c. the turn of the energy 


flows into a point is equal to the sum of the energy flows 
out of a point. Examples are: 

At point 35, 


At point 3, 


Hjq.+ ^31 ^ ^32 = ^ 


and 


H? 


H? 


At point 1: 

K? + H$j + + $4 - p i Vi = Hf + ter 

Since we assume that there is no heat transfer, no 
energy storage, and no work done in the ducts, there is no 
change in enthalpy through flow paths 5, 6, and 7. 
Consequently, enthalpy flow rate between points 5 and 12 
is designated as H§. Similarly flow path 1, flow paths 16 
and 12, path 23, and paths 25 and 31 represent constant 
enthalpy p roc esses. 

CONTROL INPUTS - The variables which tit in- 
dependent of the system and are available to be arbitrarily 
determined as a function of time or as a function of the 
system variables are designated control variables. They 
are chamber lighting, chamber wall temperatures, valve 
orifice flow areas, gas flow rates out, gas flow rates in, and 
blower motor torque. 

DISTURBANCES - Disturbances may take tbe form 
of variables which are independent of the system internal 
variables but may not be arbitrarily determined (such ms 
chamber volume changes due to glove port use) or changes 
in system parameters due to aging, etc. 

SYSTEM EQUATIONS - The constitutive relation- 
ships (processes) and structural equations form the set of 
algebraic and dynamic primitive equations which together 
describe the linked thermal and fluid dynamic behavior of 
the plant growth chamber and the atmospheric control 
system of tbe Crop Growth Research Chamber. Pressures 
and temperatures can be related to the potentials at each 
point on the molar component and energy component 
diagrams, respectively. The molar and energy flow e nd 
storage process representations shown in the proceeding 
are symbolic in that in the actual representation 
accounting for molar flow and storage by species is 
necessary. The equations can be written in terras of the 
partial pressures of each atmospheric component, total 
pressures, mass or molar flow rates (hence velocities, 
transpiration rates, etc.), air temperatures, plant canopy 
temperature, etc. 

The equations, in state variable form, can be used to 
analyse system properties such as: (I) location of 
equilibrium points, (2) stability at equilibrium points, (3) 
stability robustness at stable equilibrium points, (4) 
controllability, and (5) observability. The state variable 
form of tbe equations, formatted to include system 
uncertainties and disturbances, can be used to seek robust 
control algorithms. If analysis of the system model 
indicates that the properties of ths ,“* 

unsatisfactory such that robustness cannot be achieved or 
can only be achieved at unacceptable cost of some other 
performance measure, a proposed alteration to the system 
itrueture may be required. The model of a proposed 
reconfiguration can developed by reformulating the 

structural relationships of the process models to conform 
to structure implied by the altered configuration. 
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Fig. 24. Molar component of the modeling diagram of the Plant Growth Chamber and atmospheric 
control system of the Crop Growth Research Chamber. Points and flow paths correspond to those 
indicated on Fig. 2. Modeling symbols are defined in Figs. 3 through 23. 


Additional process models can be added if required by the over ■* growth cycle. 

altered configuration. li u characteristic of the system control problem to 

employ simplified models for the purpose of controller 
CONCLUDING REMARKS lyetheau. These model* typically (a) employ functionally 

mot* nimpic (linear, usually) representations of process 
Wn have presented a symbology for an organised behavior and (b) do not describe all stable dynamics, 

approach to tbe derivation of the equations which describe Historically, thin strategy has been remarkably tuccaatfnl 

the dynamic behavior of a plant growth chamber with because implementation of feedback control can reader a 

cloture. We have illustrated the linkage between the syntem robust to a variety of perturbation*. Computer 

model of the physical system components aad the model of simulations, based on more comprehensive and nearly 

the biological component. As previously noted, several complete models which at least include the moet 

assumptions were made to reduce the complexity of the significant nonlinearities, are most useful to demonstrate 

system for the purpose of demonstrating the modeling controlled system performance. Simulations and the 

procedure. Additional processes, such as those affecting newly 'merging perspectives of robust control are crucial 

crop shoot-root interactions, would be added in order to for the purpose of assuring performance, 

model tbe behavior of the *y*tem in response to root tone 
environment disturbances and control inputs. Proc esses 
which affect longer term phenomena, such as biomaas 
production, would be added to model system behavior 
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NOMENCLATURE 



A« 

cuopy ina (m 3 ) 

D 

duct diameter (m) 

A* 

effective croaa sectional am of tba ductwork (m 3 ) 

*« 

rate of energy accumulation at point i (waits) 

A, 

duct croaa sectional irta in flow path i (m 3 ) 

F u 

radiative shape factor between surfaces i and j 

A i 

ana of surface i (m 3 ) 

r 

friction factor 

A# 

beat transfer ana (m 3 ) 

(» 

photosynthetic rate (watts) 


total effective canopy ana for short wave 

f,i 

PAR absorption rate (watts) 


absorption (m” 3 ) 

H/ 

total enthalpy flow rate in path i (watts) 

air 

dry air 

h< 

canopy convective heat transfer coefficient 

c j 

mole fraction of carbon dioxide at point j 


(watts m” 3 K“ l ) 

** 

molar apecific heat in path i (Joules mo! -1 K“ l ) 

h j 

specific enthalpy at point j (Joules mol- 1 ) 

CO, 

carbon dioxide 

hjo 

water - liquid or vapor (in context) 
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I moment of inertia of motor end impeller (kg m 2 ) 

1^^ incident PAR above canopy {fimok m -3 sec l ) 

KA« effective area, for diffusion in the canopy (m 3 ) 


K« 


K, 

K ^j 


L 

MW. ir 


MW, 

"j 


P 1 V 1 

P i 

P j’ 

PPF 

Qi 


R 

FU 


average energy content of photoeynthetic photons 
(J jiraor 1 ) 

pump pressure reaction coefficient 
radiative heat transfer coefficient (m~ 3 ) 

K rij = K -ij ( F ij' *i* ‘j* A i' A j ) 

effective duct length (m) 

i molecular weight of dry air in chamber (g mol L ) 
molecular weight of air in the ductwork (g mol *) 
molar flow rate in path i (moles sec l ) 
nitrogen 
oxygen 

work done by the system (watts) 

pressure at point j (Pa) 

total pressure at point j (Pa) 

photosynthetic photon flux (pmoi m" 3 sec "*) 

heat transfer rate in path i (watts) 

humidity ratio at point j 

universal gas constant (Joules mol* 1 K” 1 ) 

internal fluid resistance (n m sec* 1 ) 

water vapor boundary layer resistance (sec m _l ) 
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ABSTRACT 

In support of development of bioregenerative life support 
systems for space habitats, a chamber for conducting re- 
search on the biophysical response of crop plants in con- 
trolled environments is being designed by the National 
Aeronautics and Space Administration. The imprecision 
of mathematical descriptions of the behavior of biological 
systems led to the development of a model which can oe 
used to derive a strategy for control of the chamber envi- 
ronment which is robust to the system uncertainties. The 
modeling approach and observations of the characteristic* 
of the model are described. 

INTRODUCTION 

Regenerative life support systems ars judged to be essen- 
tial to success of mission scenarios contemplated by the 
Office of Exploration of the National Aeronautics and 
Space Administration (NASA). A higher degree of 
material closure than that provided by phystochemtcai life 
support technology is considered necessary for some 
scenarios. (1] A means of providing greater closure is 
through use of a bioregenerative life support system, i.e. 
one in which the production of food, the replenishment of 
oxygen and removal of carbon dioxide from the 
atmosphere and the production of fresh water is provided 
by higher plants and/or algae in conjunction with 
physiochemicai systems. 

To conduct experiments on the biophysical res posses of 
higher plants from seed to maturation tn closed environ- 
ments. a crop growth research chamber (CGRC) [2J is 
being developed at the NASA Ames Research Center 
under the Controlled Ecological Life Support System 
(CELSS) program. [3] Stringent control of all environ- 
mental variables (e.g. lighting, temperature, atmospheric 
oxygen, carbon dioxide concentrations, humidity snd 
pressure, nutrient solution temperature, ion and dissolved 
oxygen concentrations and pressure) which affect the 
production of biomass, oxygen, and transpired water 
vapor, and the consumption of carbon dioxide, liquid 
water, and nutrients by the plants will be required. 
Knowledge of the autecology of crop plant populations 
derived from CGRC research will permit ulilisatioe of the 
CELSS plant growth unit as a multiple input-multiple 
output life support system component whose function and 
performance can be varied as required through application 
of appropriate control strategies. 

SYSTEM DESCRIPTION 

The CGRC difTers from other plant growth chambers 
(references in (4]) in size, degree of material closure, 
performance requirements, and number of controlled 
environmental variables (Table 1). Uncertainties in the 
CGRC system are moat manifest in the living component. 
The inherent variability of biological systems makes the 
task of precise mathematical description of their behavior 
difficult. The bounds set on the required performance of 


the system, coupled with the uncertainties arising from the 
imprecision in knowledge of functional forms and 
parameters describing plant biophysical ecology, have led 
us to seek a strategy for control of the chamber 
environment which is robust to the system uncertainties 
using Lynpnnov theory based methods (5). In this paper 
we describe the development of a modeling structure 
consistent with the derivation of such controllers. 

MODELING APPROACH 

The model includes the linkages between the physical 
system components and the biological components. 
Tentative physical component sizes were based upoe 
lighting, pressure, phase and temperature rise constraints 
and the proposed physical configuration of the chamber 
(Fig. 1). Assuming that the chamber contained a full 
canopy of a nearly mature crop of wheat (Yecora Rojo) 
under fully lighted conditions, essentially ail the photo- 
•yathetically active radiant energy from the lights was 
considered absorbed by the plants. A large fraction (5/6) 
was assumed converted into sensible heat, raising the ca- 
nopy temperature slightly above that of the nominal 
ambient chamber air, the remainder used in evaporating 
water internal to the plant canopy and subsequently 
transpired. The fraction of energy fixed in photosynthesis 
or released in metabolism was considered sufficiently small 
to neglect in the energy balance for short term dynamics. 
The plant canopy was modeled as temperature and vapor 
and gas exchange disturbances. A description of the 
physical system and tbe modeling system used to derive 
the constitutive and structural relationships describing the 
processes occurring in the shoot cone of tbe plant growth 
chamber is developed in (4j. 

The primitive system equations for the nominal model (i.e. 
without uncertain ties) 

4, i e A, n + # + A, « + 4, 

0| 9 4- 0} n + 0 S If + 04 4 = 0 (1) 

were derived by linearizing the process equations about an 
equilibrium point calculated by setting conditions in the 
chamber at values within the required operating envelope. 
The nominal linear time invariant slate variable form 
i * + C« + fd 

y wC*+0*i+«rf ( 2 ) 

was derived from (I) according to (6). The 19 natural 
variable system states consul of the molar compositions 
and temperatures of the air in the upper and lower 
sections of the chamber and in the ducts, blower motor 
speed, molar flow rate in the ducts and upper chamber 
volume. The LI control variables are described in terms of 
valve loss coefficients, blower motor control torque, and 
gas injection and extraction rates. The 6 disturbance 
variables are plant canopy temperature and gas exchange 
rates. chamber wall temperatures. gas injection 
temperature and chamber volume changes due to glove 
port utilization. 
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OBSERVATIONS 

Some initial obwnrttiou of th« properties of tbe nominal 
liDMr model rekvut to control design ere discussed 
below. 

a. Utilisation of the SVD in MATLAB indicated that, for 
the scenario mod el ed, B t is poorly conditioned. The 
values of some eigenvalues of fan quite sensitive to the 
method used for matrix inversion in the equation 
reduction procedure. Equation reduction from primitive 
to state variable form was performed using MATLAB. 

b. Of the 10 eigenvalues, 3 are positive real, 5 are negative 
real, one is sera, the others are complex with negative real 
parts. All pairs except one have damping ratios (O 
greater than 0.7. 

c. The s y s tem was transformed to modal form. 
Examination of the real Jordan matrix inverse revealed 
that most modes are primarily influenced by rational 
groups of natural variables. For example, the molar 
quantity of water vapor throughout the system (i.e. in the 
upper chamber, lower chamber, and in the ducts) 
predominantly influences the most unstable mode and one 
complex mode; chamber volume predominantly influences 
the second most unstable mode and a second complex 
mode; molar quantity of car boa dioxide throughout the 
system predominantly influences the ieest unstable mode 
and a third complex mode. The magnitudes of the three 
complex modes differ by less than 1%. 

d. All inodes other than that nse ori at ed with chamber 
volume are controllable. All modes are vulnerable to tbe 
disturbances. 

e. Placement of tbe largest stable eigenvalue (dominated 
by system oxygen concentration) from tbe open loop value 
of -0.257E-6 to -141 was achieved by linear feedback of 
only tbe state var ia bles representing oxygen and system 
volume without affecting more than the fifth significant 
figure of any of the other eigenvalues. 

CONTINUING DEVELOPMENT 

A paramet eri ne d linear model is being developed to 
accommodate potential variations in proponed CGRC 
design configuration. Bounds on the values for the 
uncertainties associated with the crop dynamics will be 


TABLE 1. Operating ranges and tolerance require* 
menu for the atmospheric variables in the shoot sone 
of the plant growth chamber of the NASA CELSS 
Crop Growth Research Chamber. 


Air t em p er ature 
Relative humidity 
Carbon dioxide 
Oxygen 
Nitrogen 
Gage pressure 
Air velocity 
Photosyuthetic 
photon flux 
Surface temperatures 


5-40* C ± 1* C 
35-90% ± 2%* 
25*50000 ppm ± 5%* I 
5-25% ± 5%* 
75*95% ± 5%* 
47 mm Hg ± 5 mm Hg 
0.1-1 m sec* 1 ± .1 m sec* 1 
0-3000 pmoies m _a s* 1 
± 10 pmoles m* 3 s“ l 
Air temperature + 2* C 


• tolerance ms % of set point value 


obtained from CELSS crop physiologists. Bounds on the 
values for the uncertainties associated with the physical 
lyitsn will be obtained from the NASA Ames Research 
Center CGRC engineering team. Computer simulations of 
the nonlinear system equations are being developed. The 
simulations will be used to conduct scenario studies to 
evaluate system performance with various candidate 
control strategies. 

ACKNOWLEDGEM ENT 

'Phis research was supported by NASA-ASEE Summer 
Faculty Fellowships at Stanford University - NASA Ames 
Research Center awarded to ALB and CCB- 

REFERENCES 

111 D. A. O’Handley, “Human pathway to the solar 
system in tbe 21st century,* SAE Technical Paper 
891430. 191* /ntersooety Conference on Environments/ 

Systems, San Diego, CA, 1989. 

[21 D. L. Bubenbfiim , P. M. Luna, M. Hnslerud, K. M. 
Waxen bach, and C. L. Straight, “The Crop Growth 
Rmearch Chamber, a ground bused facility for CELSS 
riTT — w " SAE Technical Paper 891588. 19t* 

/nterssciety Conference on Environmental Systems, San 
Diego, CA, 1989. 

(31 R. D. Mac Elroy and D. T. Smernoff, eds., “Controlled 
fife support system: regenerative life-support 
systems in space,* NASA CP 2480, 1987. 

(4] A. L. Blackwell and C. C. Blackwell, ‘A modeling 

system for control of the thermal and fluid dynamics of 
the NASA CELSS Crop Growth Research Chamber," 
SAE Technical Paper 891570. 19t* Mleroociety 

Con/erence on Environmental Systems, San Diego, CA, 
1989. 

(5] M Cortess **** G. Lehmann, “Controller design for 
uacertuB tyium. via Lyapunov function.." 

of ike Automatic Control Conference, pp. 2019-2024, 
Atlanta, GA, 1988. 

[81 C. C. Blackwell, “On generating state space equations 
of a linear constant coefficient system," Jovrnel of 
Opttmtxelion Theory end Affkeeitone, pp. 65-79, 1986. 



2114 




APPENDIX A. 3: 


THE LINEAR STATE VARIABLE REPRESENTATION COEFFICIENT 
MATRICES 


A3.1 



APL = 

1 . 0D+05 * 

Columns 1 thru 4 
- 0.002424541202636 
0.000007428647359 
0.000002400532336 
0.000000012012662 
0.000000158817643 
0.000000886999249 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
- 0.000007473646694 
- 0.000002399467664 
- 0.000000012002662 
- 0.000000161035711 
- 0.000071705926065 
0.000000000000000 
0.000000000000000 
Columns 5 thru 8 
- 0.015584615000000 
- 0.104155231885820 
- 0.033657271994475 
- 0.000168286354972 
- 0.002270129825046 
- 0.047268854535376 
0.104155231885820 
0.033657271994475 
0.000168286354972 
0.002270129825046 
- 0.139011870881306 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 9 thru 12 
0.000000000000000 
0.306365761428881 
0.099000649983750 
0.000495003249919 
0.006641293603077 
0.139026042751107 
- 0.306338921428881 
- 0.098991976983750 
- 0.000531103249919 


- 0.015584615000000 

- 0.104167521885820 

- 0.033657271994475 

- 0.000168286354972 

- 0.002257842025046 

- 0.047268854535376 

0.104167521885820 

0.033657271994475 

0.000168286354972 

0.002257842025046 

- 0.139011870881306 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

- 0.004425979500000 

- 0.029582335980663 

- 0.009559392265193 

- 0.000047796961326 

- 0.000641275897790 

- 0.013465337587892 

0.029582335980663 

0.009559392265193 

0.000047796961326 

0.000641275897790 

- 0.039428269392866 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.306365761428881 

0.099000649983750 

0.000495003249919 

0.006641293603077 

0.139026042751107 

- 0.306338921428881 

- 0.098991956983750 

- 0.000494960009919 


- 0.015584615000000 

- 0.104155231885820 

- 0.033669561994475 

- 0.000168286354972 

- 0.002257842025046 

- 0.047268854535376 

0.104155231885820 

0.033669561994475 

0.000168286354972 

0.002257842025046 

- 0.139011870881306 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.306365761428881 

0.099000649983750 

0.000495003249919 

0.006641293603077 

0.139026042751107 

- 0.306375061428881 

- 0.098991976983750 

- 0.000494960009919 

- 0.006640711803077 

0.408858443768547 

0.045837102941176 

0.000009300000000 

- 0.000008673000000 

- 0.000000043360000 

- 0.000000581800000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.029593147410898 

0.009562885927852 

0.000047814429639 

0.000641510264327 

0.013429105647768 

- 0.029593147410898 

- 0.009562885927852 

- 0.000047814429639 


- 0.015584615000000 

- 0.104155231885820 

- 0.033657271994475 

- 0.000180581104972 

- 0.002257842025046 

- 0.047268854535376 

0.104155231885820 

0.033657271994475 

0.000180581104972 

0.002257842025046 

- 0.139011870881306 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.306365761428881 

0.099000649983750 

0.000495003249919 

0.006641293603077 

0.139026042751107 

- 0.306338921428881 

- 0.099028119983750 

- 0.000494960009919 

- 0.006640711803077 

0.408858443768547 

0.045837102941176 

- 0.000026800000000 

0.000027470000000 

- 0.000000043360000 

- 0.000000581800000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

- 0.000007427000000 

- 0.000002400000000 

- 0.000000012000000 
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APL ( continued ) 
- 0.006640711803077 
0.408858443768547 
0.045837102941176 
- 0.000026840000000 
- 0.000008673000000 
0.000036100000000 
- 0.000000581800000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 13 thru 16 
0.031773295283869 
0.000006330000000 
- 0.000005901000000 
- 0.000000029530000 
- 0.000000395800000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
- 0.031773295283869 
- 0.000006220000000 
0.000005901000000 
0.000000029510000 
0.000000395800000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 17 thru 19 
0.004429819765695 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000012636792617 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
- 0.004429819765695 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
- 0.000024575432811 
0.000000000000000 
0.000000000000000 


- 0.006676853603077 

0.408858443768547 

0.045837102941176 

- 0.000026840000000 

- 0.000008673000000 

- 0.000000043360000 

0.000035560000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.031773295283869 

- 0.000018260000000 

0.000018690000000 

- 0.000000029530000 

- 0.000000395800000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

- 0.031773295283869 

0.000018370000000 

- 0.000018690000000 

0.000000029510000 

0.000000395800000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000519274000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000822452196 

- 0.000001580805758 

0.000000000000000 


- 0.000641510264327 

0.039441678349110 

0.004427597058824 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000024575432811 

0.000000000000000 

0.000000000000000 

0.031773295283869 

- 0.000018260000000 

- 0.000005901000000 

0.000024560000000 

- 0.000000395800000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

- 0.031773295283869 

0.000018370000000 

0.000005901000000 

- 0.000024550000000 

0.000000395800000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.648836700000000 

4.336690953038673 

1.401381215469613 

0.007006906077348 

0.094009323204420 

1.967951572080137 

- 4.336690953038673 

- 1.401381215469613 

- 0.007006906077348 

- 0.094009323204420 

5.787502839399791 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


- 0.000000161000000 

- 0.000009894665308 

- 0.000489255870000 

0.000007427000000 

0.000002400000000 

0.000000012000000 

0.000000161000000 

0.000071561445609 

- 0.000003084020638 

0.000000000000000 

0.031773295283869 

- 0.000018260000000 

- 0.000005901000000 

- 0.000000029530000 

0.000025090000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

- 0.031773295283869 

0.000018370000000 

0.000005901000000 

0.000000029510000 

- 0.000024180000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 
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BPL = 

1 . 0D+04 * 

Columns 1 thru 4 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
- 1.339113000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
Columns 5 thru £ 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000100000000000 
0.000000000000000 
0.000000000000000 
0.000725474031327 
0.000000000000000 
0.000000000000000 
Columns 9 thru 1 
- 0.000011567184703 
- 0.000000004442601 
- 0.000000001435605 
- 0.000000000007178 
0.000000005885384 
0.000000142327155 
0.000000000000000 
0.000000000000000 
0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

■ 0.000074270000000 

- 0.000024000000000 

- 0.000000120000000 

- 0.000001610000000 

- 0.000098946653075 

- 0.004591610000000 

0.000074270000000 

- 0.000076000000000 

0.000000120000000 

0.000001610000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

! 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000100000000000 

0.000000000000000 

0.000911788953009 

0.000000000000000 

0.000000000000000 

1 

- 0.533985494195872 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.001232788811738 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000017250000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000100000000000 

0.000830173124485 

0.000000000000000 

0.000000000000000 

- 0.024352741122039 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

- 0.000138158814227 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000100000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000717889530091 

0.000000000000000 

0.000000000000000 

- 0.105145008041067 

0.000000850426248 

0.000000274811229 

0.000000001374056 

- 0.000001126611533 

- 0.000027245020290 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

- 0.000000850082734 

0.000000274811229 

- 0.000000001374056 

- 0.000000018435253 

- 0.000008208843239 

0.000000000000000 

0.000000000000000 
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BPL(continued) 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000004440806 

-0.000000001435605 

0.000000000007178 

0.000000000096305 

0.000000042882747 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

1.000000000000000 

-0.428711727325246 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.101936799184506 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


CPL = 

Columns 1 thru 
0.000000000000000 
0.000000000000000 
-1.000000000000000 
0.000000000000000 
0.000000000000000 
-0.374812247255921 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 

Columns 5 thru ! 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.034662045060659 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

-1.000000000000000 

0.000000000000000 

-0.428076256499133 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

-0.077007510109763 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

2.357019064124783 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000000000000000 

0.000070873866447 
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CPL(continued) 

0.000000000000000 

0.000000000000000 


0.000000000000000 

1.000000000000000 


0.000000000000000 

0.000000000000000 


0.000000000000000 

0.000000000000000 
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